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ABSTRACT 

The  disturbance  field  generated  by  a  point  source  is 
obtained  for  steady,  three  dimensional,  isentropic,  magneto- 
hydrodynamic  flow.   The  flow  is  governed  by  the  Lundquist 
equations  linearized  about  constant  velocity,  matter  density, 
and  magnetic  field.   The  flow  is  assumed  to  be  hyperliptic, 
i.e.,  the  negative  of  the  free  stream  velocity  lies  within 
the  fast  wave  front  of  the  time-dependent  characteristic 
surface  but  outside  the  two  cusped  slow  wave  fronts  of  this 
surface.   There  are  both  real  and  complex  characteristics. 

New  dependent  variables  are  introduced  which  separate 
the  flow  into  transverse  and  compressive  components.   The 
transverse  component  is  propagated  one-dimensionally  with 
the  Alfven  velocity.   The  determination  of  the  compressive 
component  is  reduced  to  finding  the  fundamental  solution  of  a 
fourth  order  homogeneous  differential  equation  with  constant 
coefficients.   The  solution  is  required  to  satisfy  smoothness 
and  causality  conditions  and  a  condition  on  the  growth  at 
infinity. 

The  fundamental  solution  is  found  by  extending  the 
plane  wave  representation  for  the  fundamental  solution  of 
hyperbolic  equations  to  the  hyperliptic  case.   The  fundamental 
solution  for  a  hyperbolic  equation  is  extended  to  include  plane 
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waves  with  complex  wave  speeds.   This  function  will  not  be 
continuous  across  the  plane  z  =  0.   A  solution  for  z  :/  0  of 

the  homogeneous  equation  is  then  added  so  that  the  resulting 

2    2 
function  is  smooth  across  z  =  0,  x  +  y  /  0.   This  function 

is  shown  to  he  the  fundamental  solution  for  the  compressive 

component  of  the  steady  flow. 

Since  there  are  two  directions  for  which  the  real  wave 

speeds  coincide^  it  is  possible  that  the  solution  is  singular 

across  the  part  of  the  plane  ruled  by  the  tangent  half  lines 

to  the  line  segment  joining  the  slow  wave  fronts  of  the 

time-dependent  characteristic  surface.   The  steady  flow  is 

shown  to  be  smooth  across  this  surface.   The  singularities 

of  the  steady  flow  across  the  forward-facing  nappes  of  the 

characteristic  cone  are  computed. 
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Section  1.   Introduction 

We  treat  steady,  three  dimensional_,  isentropic_, 
magnetohydrodynamic  flow  past  a  point  source „   The  flow 
is  governed  by  the  Lundquist  equations  [  1  ]  which  represent 
a  non-dissipative  fluid  with  scalar  pressure  tensor _,  infinite 
conductivityj  and  in  which  displacement  currents  have  been 
ignored.   We  examine  only  small  steady  disturbances  imposed 
on  a  fluid  initially  with  constant  velocity,  matter  density _, 
and  magnetic  field;  that  is,  we  employ  the  Lundquist  equations 
linearized  about  these  quantities.   It  is  assumed  that  the 
flow  is  hyperliptic,  i.e.,  the  negative  of  the  free  stream 
velocity  lies  within  the  fast  wave  front  of  the  time- 
dependent  characteristic  surface  but  outside  the  two  cusped 
slow  wave  fronts  of  this  surface.   The  equations  for  steady 
flow  will  then  have  both  real  and  complex  characteristics. 
We  shall  find  the  fundamental  solution  for  the  steady  flow, 
i.e.,  the  disturbance  field  generated  by  a  point  source. 
The  field  quantities  must  be  regarded  not  as  ordinary  functions 
but  as  distributions.   The  fundamental  solution  will  have 
both  elliptic  properties  (e.g.  the  disturbance  will  extend 
throughout  space)  and  hyperbolic  properties  (e.g.  singulari- 
ties will  be  carried  on  the  real  branches  of  the  characteristic 
cone) . 

We  introduce  new  dependent  variables,  following 
H.  Grad  [  2  ],  which  simplify  the  analysis  by  separating 
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the  transverse  and  compressive  components  of  the  flow. 
The  transverse  component  is  propagated  one  dimensionally 
with  the  Alfven  velocity  and  its  contribution  to  the 
steady  flow  characteristic  cone  is  two  degenerate  Alfven 
cones^  i.e.  half  lines  which  are  the  loci  of  Alfven  wave 
fronts . 

The  determination  of  the  compressive  component  is 
reduced  to  finding  the  fundamental  solution  of  a  fourth 
order  homogeneous  partial  differential  equation  with 
constant  coefficients.  The  solution  is  required  to  satisfy 
three  conditions.   The  first  condition  is  a  causality 
condition  which  states  that  the  solution  should  not  be 
singular  on  the  backward-facing  nappes  of  the  character- 
istic cone.   This  type  of  condition  is  well  known  in  gas 
dynamics.   The  second  condition  is  that  the  solution  is 
required  to  be  smooth  in  all  of  space  except  on  the 
forward-facing  nappes  of  the  characteristic  cone.   Finally ;, 
we  require  the  compressive  component  to  vanish  at  infinity 

as  the  inverse  square  of  the  distance,  i.e.  as  —r^ — ?^ — py  > 

x^+y^+z^ 

except J  of  course  J  on  the  forward-facing  nappes.   We 
expect  the  fundamental  solution  satisfying  these  three 
conditions  will  be  unique;  however  we  do  not  give  a  unique- 
ness proof. 

Hyperliptic  steady  flow  has  the  property  that  the 
forward-facing  nappes  of  the  characteristic  cone  lie  in  a 
half  space,  and  we  choose  coordinates  so  that  this  is  the 


space  z  >  0.   The  causality  condition  is  then  that  the 
solution  be  smooth  for  z  <  0.   The  characteristic  polynomial 
gives  two  real  and  two  non-real  complex  conjugate  wave  speeds 
in  each  direction.   There  are  two  directions  where  the  real 
wave  speeds  coincide. 

We  find  the  fundamental  solution  by  extending  the 
plane  wave  representation  for  the  fundamental  solution  of 
hyperbolic  equations  to  the  hyperliptic  case.   If  our 
equation  were  hyperbolic,  i.e.  if  all  the  wave  speeds  were 
realj  the  fundamental  solution  would  be  H(z)u, ,  where  H  is 
the  Heaviside  function.   The  function  u-,  consists  of  the 
Laplacian  (in  x  and  y)  iterated  three  times  applied  to  a 
superposition  of  plane  wave  functions  (functions  of 
X  cos  0  +  y  sin  9   +   A.(9)z  where  'K.{9),    i  =  1,2,3^'^  is  a 
wave  speed  in  the  0-direction) .   The  plane  wave  functions 
used  include  log  \x   cos  0  +  y  sin  9   +  }\^{9)z\  ,    1  =  1,2,3,4. 
We  begin  modifying  u-,  by  extending  these  logarithms  into 
the  complex  plane  for  the  plane  waves  with  complex  wave 
speeds.   The  resulting  expression,  u„,  will  be  discontinuous 
across  the  plane  z  =  0.   We  correct  this  by  finding  a 
"smoothing"  function,  M,  where  M  is  a  solution  of  the 

homogeneous  equation  for  z  ^  0  and  u^  =  H(z)u2  +  M  is 

2   2 
smooth  across  the  plane  z  =  0,  x  +y  y'  0.   The  function  M 

is  found  by  assuming  it  to  be  a  superposition  of  plane  waves 

similar  in  form  to  u„.   However,  we  only  include  plane  waves 

with  complex  wave  speeds  so  as  to  satisfy  the  causality 
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condition^  i.e.  so  as  not  to  introduce  singularities  on 
the  backward-facing  nappes  of  the  characteristic  cone. 
The  weight  functions  included  in  the  superposition 
composing  M  are  determined  by  the  condition  that  u^  and 

its  first  three  normal  derivatives  he  continuous  across 

p   ? 
z  =  Oj  X  +y  4   0'   The  function  u^  and  its  first  two 

normal  derivatives  will  then  be  continuous  across  z  :=  0 

while  its  third  normal  derivative  1:bs  a  jump  of  a  constant 

tim.es  5(x)5(y)^  the  two-dimensional  delta  function^  across 

this  plane.   We  then  easily  exploit  this  property  to  show 

that  u,  is  a  fundamental  solution. 

The  proof  that  u^  satisfies  the  three  conditions 
above  is  straightforward.   The  causality  condition  is 
satisfied  by  virtue  of  our  method  of  selecting  M  and  we 
show  u.^  is  (real)  analytic  for  z  <  0.   We  use  the  Cauchy- 
KowalaAiski  theorem  to  show  u^  is  analytic  in  a  neighborhood 
of  z  =  0,  x^+y^  ^   0.   We  finally  show  by  a  displacement  of 
contour  method  that  u   is  analytic  for  z  >  0  except  on  the 
forward-facing  nappes  of  the  characteristic  cone.   The 
fundamental  solution  is  a  function  homogeneous  of  degree 
one.  The  compressive  component  is  composed  of  linear 
combinations  of  derivatives  of  third  and  fourth  order  of 
the  fundamental  solution  and  will  thus  have  the  desired 
rate  of  decay  at  infinity. 

The  forward-facing  nappes  of  the  characteristic 
cone  consist   of  the  two  surfaces  ruled  by  the  tangent 
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half  lines  from  the  negative  of  the  free  stream  velocity 
to  the  slow  wave  fronts  of  the  time -dependent  characteristic 
surface.   It  also  contains  the  part  of  the  plane  ruled  by 
the  tangent  half  lines  to  the  line  segment  joining  the  slow 
wave  fronts.   This  planar  segment  arises  because  of  the 
presence  of  real  double  points  in  the  characteristic 
polynomial.   The  singularity  across  it  is  computed  and  shown 
to  vanish  with  the  application  of  any  derivative  tangential 
to  the  planar  segment.  The  compressive  component  is  composed 
of  linear  combinations  of  derivatives  of  the  fundamental 
solution  and  the  fundamental  solution  is  differentiated 
tangentially  at  least  once  in  every  term  of  these  linear 
combinations.   The  compressive  component  is  thus  smooth 
across  the  planar  segment. 

For  our  method  of  solution^  we  use  coordinate  systems 
which  have  the  forward-facing  nappes  of  the  characteristic 
cone  in  the  space  z  >  0.   A  Cartesian  coordinate  system 
with  positive  z-axis  in  the  direction  of  the  free  stream_5 
such  as  is  often  used  in  steady  flow  problems,  is  then 
suitable  only  if  the  forward-facing  nappes  lie  completely 
downstream.   This  will  occur  for  equilibria  for  which  the 
negative  of  the  free  stream  velocity  lies  outside  the  slow 
normal  speed  loci.  A  Cartesian  coordinate  system  with  z-axis 
perpendicular  to  the  unperturbed  magnetic  field  has  the 
forward-facing  cone  in  z  >  0  for  every  hyperliptic  flow. 
Various  properties  of  the   characteristic  polynomial  are 
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derived  in  these  coordinate  systems. 

The  behavior  of  the  fundamental  solution  near  the 
characteristic  cone  (excluding  the  planar  segment)  has 
been  treated  in  the  hyperbolic  case^  [3  ]. 
The  behavior  here  is  similar  since  the  terms  of  the 
fundamental  solution  which  become  singular  in  each  case 
have  the  same  form.   Our  discussion  is  therefore  brief 
and  serves  chiefly  to  show  the  similarity  to  the 
hyperbolic  case. 
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Section  2.   Properties  of  the  Lundquist  Equations 

a;   The  Lundquist  Equations 

We  treat  a  perfectly  conducting,  isentropic  fluid 
described  by  a  velocity  field  u  ,    density  p,  and  scalar 
pressure  p  which  is  a  function  of  the  density  alone_,  i.e., 
P  ■=  P(p)'   The  electric  field  E  is  given  by 


r=  rx  u^ 


where  B  is  the  magnetic  field,  since  the  fluid  is  a 
perfect  conductor.   We  ignore  the  displacement  current  in 
Maxwell's  equations,  so 

\i  j^^  V  X  B ^  , 

where  J   is  the  current  per  unit  area  and  M-  is  the  specific 
inductive  capacity.   The  remaining  equation  of  Maxwell's 
equations  is 

(lA)  ^B^t  +  V  X(b'x"u)  =  0  , 

where  x\fe  have  used  E  =  B  X  u  .   Equation  (lA)  implies  that 
V«B  =  0  for  all  time  provided  it  is  zero  initially. 
We  neglect  dissipative  effects  such  as  heat 
conduction  and  viscosity.   The  equations  of  conservation 
of  mass  and  momentum  are  then 

(IB)  |^+  div  (pif)  =  0 

and 
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(IC)     P  If  +  p(u^V)u^+  a^Vp  +  ^  b'x(Vxb')  =  0 

respectively,  where  a  =  ySp/dp  is  the  speed  of  sound  and 
1 


■i-  B^X  (V  x"^  is  the  Lorentz  force  per  unit  volume, 


Equations  (l)  are  the  Lundquist  equations  for 
an  isentropic  fluid  and  form  a  nonlinear,  first  order, 
Galilean  invariant  system  for  the  functions  u,  B  and  p. 
An  equilibrium  solution  to  this  system  is  given  by 
u  =  Up|,   B  =  B-^  ,  and  p  =  p^,  where  Uq,  Bq  and  p^  are 
any  set  of  constants.   We  consider  a  subset  of  these 
constant  equilibria  which  we  will  describe  later  in  this 
section.   We  linearize  about  any  equilibrium  taken  from 
this  subset  and  look  for  steady  flows  past  a  point  source. 

B.   The  Characteristic  Cone  of  the  Lundquist  Equations 

The  subset  of  equilibria  we  consider  may  be 
described  geometrically  with  the  aid  of  the  characteristic 
cone.   This  cone  is  the  set  of  surfaces  of  singularity,  or 
wave  fronts,  which  propagate  from  a  point  source  in  the 
initial  plane  t  =  0.   We  shall  also  use  this  cone  to 
construct  the  singularity  surfaces  for  the  steady  flow  and 
to  determine  the  proper  causality  condition  needed  to  complete 
the  mathematical  formulation  of  the  steady  flow  problem.  The 
characteristic  equation  of  equations  (l)  (the  equation  for 
the  characteristic  surfaces)  will  also  be  useful.   We  give 
a  brief  summary  of  relevant  results  on  the  characteristic 
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equation  and  cone  using  [^,5,6]        and  then  describe 
the  equilibria  we  consider. 

A   characteristic  surface  <t>{x,Y ,z ,t)   =   0  for  a  first 
order  system  and  for  a  given  solution  thereof  is  a  surface 
across  which  the  solution  is  continuous  while  its  normal 
derivative  (i.e.  the  normal  derivative  of  each  component 
of  the  solution)  is  singular.   We  consider  a  point  on  the 
characteristic  surface  where  the  solution  has  the  value 
u  =  U-,  ,   B  =  B-,  and  p  -=  p-,  .   We  use  equations  (l)  to  find 
a  linear  algebraic  system  of  equations  for  the  normal 
derivative  of  the  solution  at  the  point  in  terms  of 
u -,  ,  B-,  J  and  p.  .   This  system  may  be  found  by  introducing 
a  coordinate  system  near  the  given  point  with  one  set  of 
coordinate  surfaces  given  by  (t>  (x_5y,z^  t )  =  constant^,  and 
the  other  coordinates  varying  in  each  such  surface. 
The  normal  derivative  of  the  solution  is  then  u".  ^  B,  ^  and  p^. 
This  normal  derivative  will  be  continuous  across  0  (x,y_,  z ^  t)  =  0 
unless  the  determinant  of  the  algebraic  system  vanishes. 
The  determinant  is  a  homogeneous  polynomial,  ?(*,,  4>  ,<^    ,0  ), 
in  the  components  of  the  normal  vector,  (*^_J*  ^*  ^'^    )i 
to  the  characteristic  surface  at  the  point.   The  coefficients 
of  P(<t',,(t>  ,<i>    ,<i>    )  depend  on  u -,  ,  B-,  and  p.,,  but  not  on 
the  coordinates  of  the  point,  i.e.  not  on  the  independent 
variables  (x,y,z,t),  since  the  equations  (l)  do  not  depend 
explicitly  on  the  independent  variables.   Setting  the 
determinant  to  zero,  we  get  the  characteristic  equation. 


P(<t',j0  ^  <j)  ^  ({)  )  =  0,  which  is  a  partial  differential  equation 
for  the  characteristic  surface  with  coefficients  depending 
upon  the  solution  to  equations  (l). 

If  we  linearize  equations  (l)  about  a  constant 
equilibrium,  u  =  u^,   B  =  B^,  and  p  =  p„,  and  find  the 
characteristic  equation  of  the  resulting  linear  system, 
this  equation  will  coincide  with  that  for  the  nonlinear 
system  (equations  (l))  if  we  set  u.  =  n^^    B-,  =  B^,  and  p-,=p„, 
Moreover,  the  characteristic  surface  will  be  independent  of 
the  solution  considered,  since  the  system  is  linear.   The 
characteristic  equation  for  the  linear  system  is  a  first 
order,  homogeneous,  partial  differential  equation  for  the 
characteristic  surface,  4'(xjy,z,t)  =  0,  with  constant 
coefficients  which  depend  upon  the  equilibrium  constants 


^,  B^,  and  p„.   The  characteristic  surfaces  for  linearized 


u 

steady  flow  may  be  found  by  seeking  solutions  of  the  above 


characteristic  equation  with  *  s  o. 

The  characteristic  equation  for  equations  (l), 
which  have  been  linearized  about  the  constant  equilibrium 
u  =  n^,    B  =  B^,    and  p  =  p^^,  may  be  found  by  first  writing 
equations  (l)  in  matrix  form  as 


where  the  matrices  C^,  i  =  0,1,2,3,  are  evaluated  at 
Vq  = (  B^  j  ,  and  then  changing  coordinates   as  mentioned 
before.   The  characteristic  equation  is 
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3 
dec  ,Cq*   +  ZZ  C.  *   ]  =  0 
^  ^   i=l  ^  ^i 


which,  after  a  computation,  is 


(2)  H\>^^,^y>\)    =  [*'  -  (Ao-Vo)^][(0')^-(a^+A^)(<t)'^) 


where 


(V(f))2  +  a^(A^-V*)2  (V*)^]  =  0 


'0^0 


*^  +  u^-V*,   a^  = 


•t  '  "0 
-> 


dp(pQ) 

dp 


is  the  sound  speed. 


K^ 


B 


0 


V^ 


0 


is  the  Alfven  velocity. 


Uq  =  |uq| ,  and  Aq  =  |Aq| 


The  characteristic  surfaces  for  steady  linearized 
flow  satisfy 


(3)  P(O,*^,*y,0^) 


=  [(u^-V*)^  -  (A,vV*)^] 


^0 


0 


:u".Vct,)^-(a2+A2)(u>.V4.)2(V<t>)2 


0 


+  a2(Aj-V0)2(V<D)2 


0 


The  first  component  in  equations  (2)  and  (3)  corresponds 
to  the  Alfven  or  transverse  disturbances  and  the  second  to  the 
compressive  disturbances. 

The  characteristic  cone  for  equations  (l)  linearized 
about  u  =  u„,  B  =  Bp^,  and  p  -^   p^,    at  a  given  point  is  given 
by  the  envelope  of  all  the  characteristic  planes  through  the 
point.   (Its  equation  will  be  independent  of  the  point 
considered,  since  the  characteristic  equation  has  this 
property.)   A  plane  through  the  origin  at  t  =  0  has  the  equation 
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-  Gt  +  x'^-  n^=  0  ,   x''^  (x,y,z),   n^=  {n^,n^,n^) ,    n^+ny+n^=  1, 
and  will  be  characteristic  if  (-c,n  )  is  a  root  of 

P(c,n  jH  jH  )  J   the  characteristic  polynomial. 

The  six  roots  -c.(n  ),  i  =  lj...,6j  may  be  shown  to 
be  real  for  all  n  and  equations  (l)  are  hyperbolic.   The 
envelope  of  a  family  of  characteristic  surfaces 
is  again  a  characteristic  surface,  since  the  characteristic 
equation  is  of  first  order.   The  characteristic  cone,  which 
is  the  envelope  of  the  above  planes ,  consists  of  several 
sheets  which  give  the  singularity  surfaces  for  a  point 
source  at  the  origin  at  t  ^  0.   This  cone,  or  more  specifi- 
cally, its  convex  hull,  determines  the  maximum  possible  range 
of  influence  of  the  point  source  for  the  linear  system  [  7  ]• 
The  characteristic  surface  is  defined  as  the  intersection 
of  the  characteristic  cone  with  the  plane  t  =  1.   This 
locus  is  three  dimensional  with  rotational  symmetry  about 
the  axis  of  A^^,  and  a  section  containing  "K^   is  given  in 
Figure  1  for  the  case  Aq  <  a^. 

The  characteristic  surface  when  A^-^  >  a^^  is  the  same, 
except  the  transverse  disturbances  now  lie  on  the  fast 
locus  with  A^   =  OQ  .   When  A^   =  a„,  the  two  inner,  or  slow, 
loci  join  with  the  outer,  or  fast,  locus. 

We  note  that  the  line  segment,  R' R,  joining  the  inner 
loci  must  be  regarded,  a  priori,  as  part  of  the  characteristic 
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surface^  although  it  is  not  given  by  the  envelope  construc- 
tion above.   The  reason  it  must  be  included  will  be  given 
later  when  we  discuss  the  effect  of  double  points ^  which 
are  values  of  n  for  which  two  or  more  c.(n  )  coincide. 
This  line  segment  may,  however,  be  shown  not  to  be  a 
singularity  surface  and  thus  not  part  of  the  characteristic 
surface.   ¥e  will  use  the  characteristic  cone  for 
equations  (l)  to  construct  the  characteristic  cone  for  the 
steady  flow  problem,  and  show,  by  examining  the  explicit 


solution,  that  the  line  segment  R'R  may  be  neglected  in 
this  construction. 

We  now  describe  the  subset  of  equilibria  treated. 
Consider  the  vector  -u^   from  the  origin  at  t  =  1  and  look 
at  the  section  of  the  characteristic  surface  containing 
both  -u^  and  A^.   We  denote  the  set  of  all  points  outside 
the  outer  locus  by  region  I,  those  inside  the  inner  loci 
by  region  II,  and  those  in  between  the  outer  locus  and  the 
inner  loci  by  region  III.   We  will  treat  the  subset  of 
equilibria  for  which  -u„  lies  in  region  III.   The  signifi- 
cance of  this  choice  will  be  discussed  when  we  consider  the 
characteristic  cone  for  the  steady  flow  problem  in  Section  J). 
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Section  3.   The  Steady  Flow  Equations 

A.   Derivation 

We  introduce  new  dependent  variables,  as  used 
by  H.  Grad  [  2  ];,  which  will  simplify  our  discussion  by 
separating  the  equations  into  two  independent  systems; 
one  will  have  as  its  characteristic  equation  the  first 
component  in  equation  (3)  set  to  zero,  and  the  other  will 
have  the  second  component  set  to  zero.   We  thereby  avoid 
having  to  deal  with  a  system  which  has  a  reducible 
characteristic  equation.   We  note  that  with  a  nontrivial 
geometry,  the  boundary  conditions  will,  in  general,  couple 
these  systems. 

We  linearize  the  Lundquist  equations  (equations  (l)) 
about  a  constant  equilibrium  (Bq,Uq,Pq)  and  introduce 
dimensionless  quantities  (B^u",  p)  for  the  perturbed  quantities 
(^  u";  p  ')  by  setting 


-> 


^.    B'      ->    u'         p' 

|bq|  0  Po 

The  motion  of  the  fluid  past  a  source, 

(pOPs(^'')^Po^oC('^^)^^0  ^s^^^))^   ^^  Siven  by  the  following 
equations : 

(4a)  UqV-u^  +  Uq.  Vp  =  pg 

(4b)      Ug(uQ.V)u%  a^Vp  +  AqAq  x(Vxr)  =  u^M^ 

(4C)        UqV  X(AqXu^)  +  AqIuq-V)!"^  =  AqB^ 
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where  Pg(x^)j,  M^(x^)  ,    Bg(x^)  are  C^  functions^  i.e.  they 

are  infinitely  diff erentiable  with  compact  support^  and 

B  (x  )  satisfies  V.B^(x^)  =  0.   We  also  have  the  equation 

V-B  -  0,  which  is  a  boundary  condition  at  infinity.   If 

V'B  =  0  at  infinity^  then;,  hy  taking  the  divergence  of 

equation  (^c),  we  find  V'B^^  0  everywhere. 

It  would  be  tempting;,  in  order  to  find  the  fundamental 

solution  of  equations  (4),  to  consider  the  flow  past  a  point 

source J  i.e.  to  replace  the  source  terms  p  ,  M  ,  and  B^ 

f^s  •'   s  ^      s 

in  equations  (4)  by  the  functions  p^5(x^),  M^5(x^);,  and 
B  5(x  ),    where  p^_,  M  .  and  B  are  constants.   However_,  we 
would  then  have  a  contradiction,  since,  from  equation  (4g), 

V.(uQVX(AjXu^)+AQ(u^.V)r)  =  OyAQV.r5(r)  , 


unless  B  =  0.   We  will  therefore  proceed  by  showing  that 
every  component,  v.,  of  the  solution  to  equations  (4) 
satisfies  the  equation  Lv.  =  q-,  where  L  is  the  same 
differential  operator  for  all  i  and  q.  is  a  linear  combina- 
tion of  various  derivatives  of  the  given  source  functions. 
We  then  solve  the  equation 

d°^"^5(x.) 

Lv  =.  5^(x)  6P(y)  5^(z)  ,   where  b'^lx.)  = ^-^  , 

^  dx?  ^ 

X 

and  the  solution  to  equations  (4)  will  follow  by  forming 
the  proper  superpositions.   (An  alternate  approach  for 
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resolving  the  above  difficulty  is  given  by  H.  Weitzner  [  8  ]•) 
B.   The  Compressive  Component 


We  start  by  introducing  the  new  dependent  variables 

> 

0' 


'a,^,^),    following  H.  Grad  [  2  ]>    with  a  =  V-u  ,    p  --  ^n''^ 


and  Y  =  A^-    B  .   These  three  variables  and  the  density  p 
satisfy  the  following  equations.   Equation  (4a)  gives  directly 


(5A)  UqQ  +  Uq  •  Vp  .=   p^  . 

The  divergence  of  equation  (4b)  gives 

(5B)      UQ(uQ.V)a  +  a^Ap  +  Aq  Ay  =   u^V  •  M^ 

The  inner  product  of  A"^  with  equations  (4b)  and  (4c)  gives 


(5C)      UqU^.  Vp  +  a^A^.  Vp   =   UqAJ.M^ 


and 


,5D)      u^A^a  -  UqAJ.VP  +  AquJ.VY  .  Aq^q  *  ^s 


These  equations  are  equivalent  to  a  single  fourth 
order  equation  for  each  dependent  variable  of  the  form 
Lv.  =  q.,  as  indicated  above.   These  equations  may  be  found 
by  first  writing  equations  (5)  in  matrix  form  as 
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— > 


where  A.    are  constant  matrices^  v 

}  os  \    rh\ 


a 


and 


r  = 


Rr 


R-. 


\ij 


We  may  find  the  equation  satisfied  by  each  component 
in  a  systematic  way  "by  formally  solving  the  algebraic  system 


1  X 


'2  y 


:)  z- 


where  the  differential  operators  D  ,  D  ,  D  are  treated  as 
constants  [7  ].   The  solution  is 

det  (A^D^+  A^Dy+A^D^)?^^  (A^D^  +  A2Dy  +  A-^D^  f  r" 

where  the  matrix  of  the  transposed  cofactors  of  the  elements 
of  a  matrix  B  is  denoted  by  B  .   Each  component  of  v   then 
satisfies  the  same  equation  but  with  a  different  forcing 
term,  namely'. 


4 


h 


(6)   det  (A.D^  +  A^D  +  A^Djv  =  y2_  2 


1  X 


2  y   p  z^  1 


b  .   D,    R. 
imn  Kmn   j 


operator 


where  b.   are  constants,  and  D. 


j=l  £+m+n=^ 

i  =  1,2,3,4 
is  the  differential 


In  other  words,  the  right  hand  side  of  equation  (6) 
consists  of  a  linear  combination  of  the  partial  derivatives 
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of  the  R. J  i  =  1,2,5,4,  where  the  degree  of  each  term  is 
three  or  four.   We  will  explicitly  give  some  of  the  forcing 
terms  when  we  discuss  the  singularity  surfaces  of  the  solution 
in  Section  6-B. 

It  will  then  be  sufficient  to  solve  the  equation 

(7)    det  (A^D^+  A2Dy+  A-jD^)v  =  5«(x)  b^iy)      5^(z) 

where  a  +  p  +  Y=6   or  7. 

We  then  obtain  the  solution  to  equations  (5)  by  forming  the 
proper  superpositions  of  these  solutions,  using  various 
values  of  a,  p  and  y. 

The  characteristic  equation  for  equations  (7)  (their 
characteristic  equations  are  identical  since  they  do  not 
depend  on  the  forcing  term)  is  the  same  as  that  for 
equations  (5).   (The  characteristic  equation  for  a  single 
higher  order  equation  is  defined  in  essentially  the  same 
way  as  for  a  system.)   Moreover,  it  may  be  obtained  from 
equation  (5)  by  setting  the  compressive  component 
to  zero. 

C.   The  Transverse  Component 

We  obtain  two  additional  equations  for  the  new 
dependent  variables  j  =  Aq«(VX  b'')  and  o)  =  Aq'(VX  u""), 
again  following  H.  Grad  [  2  ].   These  equations  are  found 
by  applying  Aq* curl  to  equations  (4b)  and  (4c) .   The  results 
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are,  denoting  the  directional  derivatives  A~   A^'V  and 
u~  Uq'V  by  d/dp  and  d/dq  respectively, 

and 

(8B)      "0  3f  -"oi  =  *0-<'-B^)  ■ 

Combining  these  equations,  we  find  that  both 
J  and  0)  satisfy  one -dimensional  wave  equations  with  wave 
speed  Ap|/u„  (the  Alfven  speed  A^  in  our  original 
coordinates),  namely  I 

2 
^0^2    ^0  ^  2   ~   ^0  dq  ^^0     %^         Uq  dp  ^^0     ^s  ^ 

^0^2    ^0  ,  2   ^   ^0  dp  ^^0     ''S^         ^0  dq  ^^0  ^   ^s ''  * 
dq       dp  -^  ^ 

The  characteristic  equation  for  the  system  (8)  is 
that  obtained  by  setting  the  transverse  factor  in  equation  (j) 
to  zero.   The  solution  of  the  above  wave  equations  is  well 
known,  and  we  shall  be  concerned  in  the  sequel  only  with 
equations  ( 5) . 

D.   Inverting  the  Solution 

After  the  solution  of  equations  (5)  and  (8),  there 
remains  the  problem  of  solving  for  u  and  B   given 
ct,  ^,    y,    j,    and  0)   (which  are  V-u^,  Aq*u^,  Aq'B'',  Aq'(VXb'*), 
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and  A^*(VXu^),  respectively).   Since  the  components  of 
u^   and  B^  along  Ap^  are  known,  we  have  four  differential 
equations  (three  from  inverting  a,  j,    and  ^,    and 
the  equation  div  B^=  O)  to  determine  these  functions  on 
the  plane  orthogonal  to  Aq.   H.  Grad  [  2  ]  showed  that  these 
equations  are  invertible  if  a  houndedness  condition  is 
imposed  at  infinity  and  if  the  inhomogeneous  terms, 
cc>    Pj  Yi  j;i  a^^d  <^are  sufficiently  smooth.  In  our  case, 
we  may  directly  invert  the  appropriate  linear  comhinations 
of  the  solutions  of  equation  (7)  without  the  necessity  of 
first  obtaining  smooth  source  functions,  since  these 
solutions  will  be  given  as  various  derivatives  applied 
to  smooth  functions.   We  may  first  invert  these  smooth 
functions  and  then  apply  the  derivatives.   This  procedure 
is  correct  since  derivatives  of  distributions  may  always 
be  interchanged. 

A  full  discussion  of  the  implications  of  these 
global  resolutions  (with  a,  p.,  Y,  j,  and  cd  smooth)  is 
given  by  H.  Grad  [  2  ]. 

E.   Construction  of  the  Steady  Flow  Characteristic  Cone 

We  may  use  the  characteristic  surface  for  equations 
(1)  to  give  a  geometric  construction  of  the  characteristic 
cone  for  equations  (5)   [6].   To  completely  determine  the 
solution  of  equations  (5),  it  is  necessary  to  add  a 
causality  requirement.   This  may  be  obtained  by  considering 
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the  steady  flow  to  be  the  limit   in  infinite  time^  of  a 
time  dependent  flow.   Equivalently,  since  equations  (l) 
are  Galilean  invariant,  we  may  regard  the  steady  flow  past 
a  point  source  as  the  result  of  the  source  moving  with 
velocity  -u^  from  infinity  into  a  stationary  fluid. 
At  each  point  in  time,  the  source  causes  a  wave  pattern, 
as  shown  in  Figure  1,  to  propagate.   The  possible 
singularity  surfaces  for  the  steady  flow  are  given  by  the 
envelope  of  these  waves  for  t  <  t^.      This  gives  a 
single-napped  (forward-facing)  cone  while  equations  (5) 
will  have  the  full  cone  (both  nappes)  as  its    characteristic 
cone,  i.e.  as  its  set  of  possible  singularity  surfaces.   We 
therefore  impose  the  causality  condition  that  the  solutions 
we  seek  be  smooth  across  the  backward-facing  nappes.   The 
forward-facing  cone  may  be  constructed  by  drawing  the 
vector  -u^  from  the  origin  at  t  =  1  and  then  the  forward- 
facing  cone  is  ruled  by  the  tangent  half  lines  from  -Uq 
to  the  characteristic  surface. 


There  are  three  qualitatively  different  cases 

-> 

^0 


depending  on  whether  -u^^  lies  in  region  I,  II  or  III. 


Figure  2  shows  the  plane  containing  -u^^  and  Aq,  with  -Uq 
lying  in  region  III.   For  our  case,  -u^  lies  in  region  III 
and  the  construction  gives  a  cone  (single-napped)  which 


is  composed  of  two  branches  ruled  by  the  tangent  half  lines 
from  -u^  to  the  two  cusped  slow  loci  of  the  characteristic 
surface   of  equations  (l),  as  well  as  an  additional  planar 
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segment  composed  of  tangents  to  the  line  joining  the  slow 
loci.   As  noted  above^  we  shall  show  that  this  planar  segment 
is  not  a  singularity  surface  for  the  steady  flow. 

The  characteristic  cone  may  also  be  found  by  forming 
the  envelope  of  all  characteristic  planes  through  the  point 
source,  and  then  using  the  causality  condition  to  choose 
the  proper  nappes.  We  set  n^=  (cos  9,    sin  B)    and  x  =  (x,y). 
A  plane  with  equation  x^.  n^-  Az  =  0  is  characteristic  if 
-A(0)  is  a  root  of  the  characteristic  polynomial  of 
equations  (7).   The  characteristic  polynomial  is 


;cos  0,sin  Q ,'h)    =  [Uq'(cos  0,sin  0,A)] 


(9)    -  (a^+Ag)[uQ-(cos  e.sin  6,\)f{l+\^) 


+  a^   (Aq-(cos  9, sin  0 ,\))       (l+A  ) 


Each  real  root  -a(9)  (0  ji  9  ^  2Tr)  of  this 
polynomial  gives  rise  to  a  branch  of  the  cone.   We  would 
then  expect  the  polynomial  to  have  two  real  and  two  nonreal 
roots  in  region  III.   We  shall  in  fact  be  able  to  prove  this 
algebraically  for  a  large  part  of  region  IIIj  for  two 
specified  coordinate  systems.   The  planar  segment  of  the 
cone  is  not  given  by  the  envelope  construction  but  arises 
due  to  the  presence  of  real  double  points  in  the  character- 
istic polynomial;,  i.e.  points  where  the  real  roots  coincide. 
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This  will  be  discussed  when  we  analyze  the  solution  in 
Section  6-A. 

The  method  used  to  solve  equations  (7),  and  hence 
equations  (5),  requires  the  existence  of  an  open  half  space 
containing  the  forward-facing  single-napped  cone,  and  we 
introduce  coordinates  so  that  this  is  the  space  z  >  0. 
We  can  (nonuniquely)  accomplish  this  for  any  point  in 
region  III  as  follows!   through  the  origin  of  the  steady 
flow  coordinate  system  (which  in  Fig.  2  is  at  -Uq)  we  take 
the  positive  z-axis  of  a  Cartesian  coordinate  system  as  the 
half  line  perpendicular  to  the  direction  of  Aq  and  inter- 
secting the  line  containing  Aq.   All  half  lines  tangent  to 
the  time  dependent  characteristic  surface,  and  hence  the 
single-napped  cone,  will  lie  in  z  >  0.   We  note  that  this 
property  even  holds  for  those  equilibria  in  region  III 
which  have  part  of  their  forward-facing  cone  upstream. 
This  property  also  holds  in  region  I,  but  never  in  region  II. 

The  z-axis  for  our  case  (-u^  in  region  III)  is 
somewhat  analogous  to  a  time  axis  for  a  non-steady  problem, 
in  that  the  discontinuities  of  the  solution  propagate  with 
increasing  z.   The  analogy  is  not  complete  in  that  we  expect 
the  solution  to  extend  throughout  space,  due  to  the  existence 
of  nonreal  roots  of  the  characteristic  polynomial.   The  cone 
then  does  not  give  the  range  of  influence  for  the  source. 
The  mathematical  formulation  of  our  problem  has  thus  both 
hyperbolic  and  elliptic  behavior,  and  the  flow  is   called 
hyperliptic . 
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The  case  of  -u^  lying  in  region  I  has  been  treated 
elsewhere  [  9  ].      The  forward-facing  nappes  lie  in  a 
half  space  and  (in  a  coordinate  system  with  these  nappes 
in  z  >  O)  the  characteristic  polynomial  has  all  real  roots. 
The  f oiTward-facing  nappes  are  composed  of  three  branches: 
an  outer  branchy  ruled  by  half  lines  to  the  fast  locus, 
and  a  branch  to  each  of  the  slow  loci.   Here,  as  in  our 
case,  the  solution  is  smooth  across  the  planar  segment 
ruled  by  the  tangent  half  lines  to  the  line  connecting  the 
slow  loci.   The  mathematical  formulation  of  the  problem  in 
region  I  leads  to  a  Cauchy  problem  for  a  hyperbolic  equation, 
and  the  z-axis  is  completely  analogous  to  a  time  axis.   In 
particular,  the  convex  hull  of  the  characteristic  cone  gives 
the  maximum  range  of  influence  of  the  point  source. 

For  -u^   in  regionll,  the  forward-facing  cone  consists 
of  two  branches;  however,  it  does  not  lie  completely  in  a 
half  space.   Ttie  mathematical  formulation  will  thus  differ 
from  that  in  regions  I  and  III,  and  this  problem  remains  to 
be  treated. 

F.  A   Note  on  Uniqueness  of  the  Steady  Flow 

We  close  this  section  by  noting  that  we  do  not 
prove  that  the  steady  flow  solution  obtained  is  unique, 
although  with  the  causality  condition,  as  well  as  the  growth 
at  infinity  and  smoothness  conditions  given  in  Section  5--A, 
uniqueness  seems  likely.   One  difficulty  is  that  the  nature 
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of  the  characteristic  polynomial  (equation  (9))  may  be 
qualitatively  different  in  different  coordinate  systems. 
For  example^  if  a  characteristic  plane  is  chosen  to  be  the 
plane  z  =  0^  we  find  an  infinite  root^  i.e.  the  polynomial 
will  be  of  lower  order  for  the  corresponding  value  of  9. 

We  shall  assume  that  the  causality  condition  together 
with  the  smoothness  and  growth  conditions  to  be  given 
(Section  5-A)  are  sufficient  to  insure  the  uniqueness  of 
the  solution  of  equation  (7)  and  hence  of  equations  (5). 
We  may  therefore  use  any  appropriate  coordinate  system 
to  facilitate  the  analysis.   We  consider  two  such  systems 
in  the  following  section^  and  derive  properties  of  the 
characteristic  polynomial  in  each  of  these  systems. 
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Section  4.   The  Characteristic  Polynomial 

We  have  reduced  the  problem  of  finding  the  steady 
flow  to  that  of  finding  a  solution  to  equation  (7)  which 
satisfies  various  conditions  to  be  described  in  Section  5-A . 
In  this  section,  we  derive  certain  properties  of  the 
characteristic  polynomial,,  which  will  be  used  later  to 
derive  and  analyze  the  steady  flow.   We  use  the  coordinate 
system  given  in  Section  5-E  which  has  the  forward-facing 
cone  in  z  >  0  for  all  equilibria  in  region  III.   A  second 
coordinate  system,  commonly  used  in  gas  dynamics,  is  also 
discussed.   Tliis  coordinate  system  has  the  forward-facing 
cone  in  z  >  0  for  a  subregion  (III-A)  of  region  III. 

We  conjecture  that  the  following  properties  of  the 
characteristic  polynomial  are  true  in  any  coordinate  system 
which  contains  the  forward-facing  cone  in  the  space  z  >  0. 
The  characteristic  polynomial  Q,  (cos  Q,    sin  9,  A)  (equation  (9)) 
will  have  two  real  and  two  non-real  (complex  conjugate)  roots 
for  every  real  9.      Moreover,  there  will  be  only  two  double  points, 
9  =  9^  and  9  =  9„  +  7r,  0<_9<_2f,  where  the  real  roots  coincide. 
We  expect  that  there  are  always  two  real  roots  since  the 
characteristic  cone  has  two  real  branches.   If  there  were 
more  than  two  real  roots  on  some  9  interval,  these  two 
branches  would  be  doubly  covered  for  those  values  of  9, 
while  essentially  the  same  branches  for  -u„  in  region  I 
are  singly  covered.   Finally,  we  shall  see  that  a  real  double 
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point  9    ,    0   <_  9^   <  w,    will  cause  a  planar  segment  of  the 
characteristic  cone  to  appear.   However,  our  construction 
only  gives  one  such  segment^  so  we  believe  there  is  only  one 
real  double  point. 

We  need  only  consider  the  polynomial  for  0_<9<7r,  for 


if  7\(0)  is  a  root,  we  immediately  have  -A(9)  and  -A(9)  are 
roots  at  the  point  9   +Tr.   We  choose  the  latter  to  give 


(10)  A(0+7r)   =   -  A(0)  . 

The  conjugate  is  chosen  since  the  complex  roots  never 
cross  the  real  axis.  We  remark  that  our  method  of 
solution  is  easily  carried  over  to  the  case  of  any  finite 
number  of  real  double  points. 

We  now  consider  the  coordinate  system  described  in 
Section  3-E,  and  we  additionally  specify  Aq  ^^  Ap|(0,l,0) 
and  u^  =  Up^(0,  cos  *,  sin  *).   The  characteristic 
polynomial  (equation  (9))  is  then 

Q,(cos    0,sin   9,a)=u„(cos    0   sin  9  +  A   sin   4>) 

(11) 

+  u{i(l+A2)[   ^^^-|  -    (^+i)    (cos    0   sin  9   +  A   sin   ^f] 

m  M  m        M 

where  m  =  ''^n^^O  ^^^   ^  '^'  ^O^^O  ^^'^   ^^^  Mach  numbers. 

We  first  show  that  Q,  always  has  two  real  roots. 
Figure  1  shows  that  a  necessary  condition  for  -Uq  to  lie 
in  region  III  is 
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1  .  1      1 
—  +  —   >   1  . 

ra    M 

We  will  not  consider  the  degenerate  case  of  parallel 
flow  where  sin  *  =  0  and  the  characteristic  polynomial  is 
of  second  order. 

The  coefficient  of  A   in  equation  (11)  is 

u^  sm  4)  (sm  *  -  (— p  +  — p) )  :, 

m    M 

which  is  negative^  by  inequality  (12) ^  so  Q  is  always  of 

fourth  order. 

Inequality  (12)  further  implies  that  Q  will  always 

be  negative  at  A  =  +oo,  but  since  Q,  is  nonnegative  when 

A  =,  -  cot  4>  sin  9,  where  Q(co3  Q ,    sin  Q,    -cot  *  sin  Q) 

uT^sin^e       2     P 
=   ^p  p —  (1+cot  4)  sin  9);,  Q,  always  has  two  real  roots, 

m  M 
Moreover^  this  shows  that  at  any  point  where  Q  has  two 

real  and  two  nonreal  roots ^  the  real  roots  can  only  coincide 

when  Q  (cos  Q,    sin  Q,    -   cot  *  sin  9)  =  0.   This  can  only 

occur  for  sin  9=0,  where,  as  we  shall  see  below,  there 

is  a  double  point  with  both  real  roots  vanishing. 

A  proof  that  there  are  two  complex  roots  in  region  III 

would  be  more  complicated,  since  these  roots  are  real  for 

some  equilibria  in  the  other  regions.   For  example,  they  are 

real,  in  this  coordinate  system,  provided  -u„  lies  in  the 

I  I     /'2   2 
subdomain  of  region  I  given  by  jzl  >  yB.^-\-k^      (which  implies 

I  sin  $1  >  V— p  +  — 9  )•   We  shall  instead  verify  the  existence 
V  m^    M 
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of  two  nonreal  roots  in  several  limiting  regions. 

We  begin  with  the  case  sin  0=0  for  all  points  in 
region  III^  since  this  value  of  9   will  not  he  covered  in 
our  various  limiting  regions  (because  of  the  inversion  of 
limits).   Without  loss  of  generality^  we  take  0  <_  9  <  tt 
and  0  <  *  <  27rj  o  y'  f.   The  roots  for  sin  0=0  may  be 
found  exactly,  since  the  polynomial  is  then 

Q(1,0,A)  =  A  u'Z  sin^*-  uj.(-^  +  -^)    (1  +  A^)A^  sin^  4) 


^0        -OV    m2' 


with  roots 


A^  =  0  and 

7,2  ^    1/m^  +  1/1M^      ^  Q 

(sin2$-  (4+^)) 
m"^   M 

by  inequality  (12),  so  our  condition  is  verified  here,  i.e. 

two  roots  are  real  and  two  complex.   We  also  have  exhibited 

the  doable  real  root  A  =  0  for  sin  0=0.   Moreover, 

there  must  be  two  nonreal  roots  in  a  neighborhood  of 

sin  0  =  0,  by  the  continuity  of  the  roots  in  sin  0  and  since 

the  complex  roots  must  coincide  to  become  real.   The  other 

two  roots  must  then  be  real  for  sin  0  small,  since  two  real 

roots  always  exist.   In  our  remaining  calculations,  sin  0 

will  always  be  bounded  away  from  zero. 

We  now  consider  the  case  of  perpendicular  flow,  i.e. 

cos  4)  =  0,   Here  the  roots  may  again  be  found  explicitly. 

In  this  case,  the  polynomial  is 
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Q(cos  6,    sin  Q,    a)=u^(1-(^  +  J_) )  A^ +u!t(%^-2^  -{\+   i-))^^ 

m    M         m  M     m   M 

ujsin  9 

+     2~2~ 

m  M 

2  2 

which  is  second  order  in  A  ,   It  is  negative  for  A  =  +  od 

2 
and  positive  for  A   =  0  so  vje  have  one  positive  and  one 

2 
negative  root  A  ,    which  gives  two  real  and  two  nonreal 

roots.   This  result  also  holds  in  some  neighborhood  of 

cos  *  =  0,  by  the  previous  argument. 

We  next  consider  the  region  near  parallel  flow^  i.e. 

sin  4)  =  0.   There  are  always  two  real  roots  for  sin  4)  y'  0, 

and  it  is  interesting  to  investigate  their  behavior  at 

sin  *  — >  0.   Expanding  the  roots  for  sin  *  small^  we  find 

two  real  roots  which  behave  to  leading  order  like  1/sin  *. 

This  follows  since,  for  A  =  0(l/sin  *),  the  leading  order 

term  satisfies 

2 
A^[  -— p-  -  (^  +  i)  (cos  0  sin  9  +  A  sin  0)^]  =0. 
m  M      va"^        M 

We  discard  the  zero  roots  -  the  two  remaining  roots 
may  be  trivially  found  and  are  real  and  distinct. 

The  remaining  roots  for  sin  *  small  are  of  zero  order 
in  sin  4)  and,  by  continuity,  will  behave  like  the  two  roots 
for  parallel  flow.   These,  as  is  well  known,  are  complex  in 
region  III^  so  our  result  is  true  near  parallel  flow. 

We  may  also  show  this  result  for  u„  small  (i.e.  for 
m,M  <<  l)  by  the  same  methods.   We  again  find  that  two  real 
roots  become  infinite  and,  to  leading  order,  the  two  zero 
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o       2 
order  roots  satisfy  sin  0(1 +A  )  ^  0,  which  gives  the 

result  here. 

As  noted  before^  the  only  double  point  for  an 
equilibrium  where  Q  always  has  two  real  and  two  nonreal 
roots  is  at  sin  0  =  0^  where  the  real  roots  vanish.   The 
corresponding  characteristic  plane  is  x  =  0^  which  is  the 
plane  of  u„  and  A^. 

Another  useful  coordinate  system  has  its  z-axis  in 
the  u'^   direction.   We  take  A^   =  Ap|(OjSin  4>jCOS  <t>)    and 
u^  =  u„(0^0,l).   The  angle  4>  is  then  the  angle  between 
u^  and  A^  (in  this  coordinate  system  as  well  as  in  the 
one  above).   The  characteristic  polynomial  (equation  (9)) 
in  this  coordinate  system  is 

Q(cos  e,sin  9,-k)  h         i         i        q  o 

n =  >^  -  (^+  ~)A^(1+A^) 

/   N      ut;  m"^  yf 

(13)       0 

+  -^-p  (1+A~)  (sin  0  sin  0  +  A  cos  <t) ) 

This  coordinate  system  has  the  disadvantage  that 
only  points  in  a  subregion^  region  Ill-A^  of  region  III 
will  have  the  forward-facing  nappes  of  their  characteristic 
cone  in  the  space  z  >  0.   In  the  plane  of  u^  and  Aqj  the 
boundary  of  region  III-A  is  given  by  those  equilibria  which 
have  their  y-axis  tangent  to  the  slow  loci  of  the  character- 
istic surface.   An  equilibrium  will  then  lie  on  the  boundary 
of  region  III-A  if  -Uq  is  the  closest  point  to  the  origin  on 
any  line  tangent  to  the  slow  loci.  The  boundary  is  then  the 
pedal  curve  of  the  slow  loci.   (The  above  condition  is  the 
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definition  of  a  pedal  curve.)  This  pedal  curve  is^ 
however^  just  the  slow  branches  of  the  normal  speed  locus, 
which  is  the  polar  plot  of  the  slow  compressive  wave  speeds 
in  each  direction.   This  locus  is  well  known  and  is  shown, 
along  with  the  fast  branch,  in  Figure  5^  where  it  is 
superposed  on  the  characteristic  surface.   In  fact,  the 
characteristic  surface  is  usually  found  as  the  envelope  of 
all  lines  which  intersect  the  normal  speed  locus  at  their 
closest  distance  to  the  origin.   This  is  the  inverse  of  the 
construction  which  gives  the  normal  speed  locus  as  the 
pedal  curve  of  the  characteristic  loci.   The  boundary  may 
also  be  found  analytically,  and  consists  of  those  equilibria 
for  which  the  characteristic  polynomial  is  of  lower  order,  i.e. 

(14)         (1  -  (%  +  4-)  +^)  =     0  . 

m    M     m  M 

The  region  III-A  now  consists  of  those  points  in 
region  III  which  lie  outside  the  two  branches  of  the  slow 
normal  speed  locus.   These  points  are  precisely  those  whose 
single-napped  cone  lies  completely  downstream.   (Downstream 
in  this  coordinate  system  is  just  the  space  z  >  0,  and 
upstream  is  z  <  0. )   Points  in  the  remainder  of  region  III 
(and  those  points  in  region  I  which  lie  within  the  fast 
branch  of  the  normal  speed  locus)  will  have  their  single- 
napped  cones  lying  both  upstream  and  downstream.   If,  and 
only  if,  the  normal  speed  loci  and  characteristic  loci 
coincide  (as  for  the  wave  equation)  will  all  equilibria  have 
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their  forward-facing  cones  downstream. 

The  proof  that  Q,  (equation  (13))  has  two  real  roots 
is  similar  to  the  proof  given  earlier.   Equation  (l4)  gives 
the  boundary  of  region  III-A  as  the  set  of  all  points  for 
which  the  highest  coefficient  of  Q  vanishes.   We  then  see 
"by  checking  a  point  inside  region  III-A,  such  as  a 
pei*pendicular  flow,  that  for  points  in  region  III-A,  Q  has 

its  leading  coefficient  negative.   Q  is  nonnegative  for  A=  0, 

u7,sin  *  sin  0 
where  it  equals  —^ p-p ,  so  Q,  has  at  least  two  real  roots 

m^M 
This  implies,  by  the  same  argument  as  above,  that  the  only 

real  double  point  must  be  at  sin  9=0,  if  Q  has  two  nonreal 


roots  for  all  9.      The  characteristic  plane  corresponding  to 

Tq  and  Aq, 


this  double  point  is  x  =  0.   This  is  the  plane  of  u^^  and  A^, 


as  we  found  in  the  other  coordinate  system. 

The  polynomial  has  two  nonreal  roots  for  perpendicular 
flow,  i.e.  cos  *  =  0.   In  this  case,  the  polynomial  is  second 
order  in  A  and,  by  exactly  the  same  argument  we  used  for 
perpendicular  flow  in  the  other  coordinate  system,  we  find 
two  real  and  two  nonreal  roots.  The  result   is  then  true  in 
a  neighborhood  of  cos  ^  =  0,  by  the  previous  argument. 

We  now  show,  following  Rellich  [10],  that  the  roots 
of  the  characteristic  polynomial  are  analytic  in  9   for 
0  <  9  <  27r   if  the  polynomial  has  two  real  and  two  nonreal 
roots  for  all  9   and  only  a  finite  number  of  real  double 
points.   The  only  difficulty  is  at  a  point  9^   where  the 
real  roots  coincide,  since  the  roots  of  a  polynomial  with 
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analytic  coefficients,   noncoincident  roots,   and  leading 
coefficient  nonzero  are  analytic.       Let  the  real  roots 

be  A^  ^{e)    with  ^^(9^)  =  ^^(^o)'  '^^^   functions  A^  ^{9) 

1/2 
have  expansions  in  {0-9    )    '      about  9^  since  the  complex 

roots  are  bounded  away  from  the  real  axis.       The  terms 

n  /2 
(0-9^)  '    ,    n  odd,  do  not  appear  in  these  expansions,  since 

A-]_  ^{9)    are  real  for  both  0-0q  >  0  and  9-Q^   <_   0.   For  if 

A^(e)  =  ;> ®n/2^^"^0^   ^  ^^^   ^m/2>    ^   °^^^  ^^  ^"^^  first 

nonzero  coefficient  of  the  form  a  /p,  n  odd,  then 

^~^  n  /P 

h(«)  -^=n/2(e-''o)  "^^ 

n=0    ' 

is  a  real  function  for  9   real.   However,  its  imaginary  part, 

to  leading  order  for  (9-9  )  small,  will  be  Im  ^r^/o^^~^ q)'"^      > 

which  will  be  nonzero  for  either  9-9      >  0  or  9-9      <   0, 

unless  a  /„  =  0,  proving  the  result. 
m/2     •'  '^  ^ 

We  close  this  section  with  a  useful  simple  property 
of  the  characteristic  polynomial  in  any  Cartesian  coordinate 
system.   By  equation  (lO)  and  the  explicit  form  of  Q, 
equation  9,    we  find 

Q,^(cos  (9+7r),  sin  {9+Tr) ,    A^(0-fir)) 


(15)  =   Q,(-  cos  9,  -  sin  9,  -  A^(9)) 


-  Q^(cos  9,  sin  9,  A^(9))  , 

1  =  1,2,3,4. 


-34- 


Section  5.   The  Derivation  of  the  Solution 
A.      Statement  of  the  Problem 

In  this  sectionj  we  find  a  solution  to  equation  (7) 
which  satisfies  the  causality  condition  and  the  smoothness 
and  growth  at  infinity  conditions  given  below.  The  steady 
flow  past  a  source  is  then  given  by  forming  the  proper 
superpositions  of  this  solution,  using  various  values  of 
a,  ^,  and  y,  and  proceeding  as  shown  earlier.  We  begin  by 
describing  the  differential  equation  and  the  conditions  we 
impose  on  its  solution. 

The  differential  equation  has  the  form 

(16)  q{b/dx,    ^/dy,    a^z)u(x,y,z)  =  5«(x)  6p(y)  5^(z)  , 

a+p  +  Y=5  or  7,  a^p.Y^l^ 

where  Q,  is  a  homogeneous  polynomial  of  fourth  degree  with 
constant  coefficients  which  depend  upon  the  Cartesian 
coordinate  system  chosen.   Such  properties  of  equation  (I6) 
as  the  number  of  roots  A(0)  of  its  characteristic  polynomial 
Q(cos  9,    sin  9,    A)   (i.e.  whether  or  not  Q(0,0,l)  vanishes) 
and  the  number  of  real  roots  of  this  polynomial  are  not 
invariant  under  arbitrary  rotations  of  the  coordinate  system. 
We  have  shown  that  in  a  large  part  of  region  III  (and  possibly 
in  all  of  region  III),  we  may  choose  at  least  one  coordinate 
system  for  each  point,  i.e.  for  each  equilibrium,  such  that 
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equation  (l6)  has  the  following  properties!  Sq  =  Q(0, 0,1)^0, 
and  for  every  unit  vector  n  =  (cos  6,    sin  9),    the  character- 
istic polynomial  Q(cos  9,    sin  9,    A)  {=   Q(n^?s))   has  roots 
-A^(n'')   (5  -A^(e)),  i  =  l,2,^,h      with  A^^2^^^)  ^^^^  ^"^^ 
A^  ii(n'')  complex,  and  Im  A^  >  0.   The  roots  all  satisfy 
A^(9  +  7r)  =  -Array  and  Q^  (ET^^-fn^), -A^  (9+7r) )  =  -Q^(HTeT. -^i(9) )  • 
The  roots  are  moreover  distinct  except  at  two  double  points, 
9   =   9      and  9   =   9^+it   ,  0  <_  9  <  2tt ,    where  the  real  roots 
coincide.   In  addition,  the  coordinate  system  can  be  chosen 
so  as  to  place  the  single-napped  forward-facing  character- 
istic cone  given  by  the  construction  of  Section  3-E  in  the 
space  z  >  0.   The  causality  condition  then  takes  the  form 
that  the  solution  be  smooth  for  z  <  0.   For  def initeness, 
we  will  assume  that  the  coordinate  system  chosen  is  one  of 
the  two  explicitly  given  in  Section  4,  which  places  the 
double  point  at  sin  9  =  0  where  A-^  ^   =   0.   We  note  that  the 
arguments  we  shall  use  may  be  extended  simply  to  the  case 
of  a  finite  number  of  double  points  where  the  real  roots 
coincide . 

We  require  the  solution  of  equation  (l6)  to  be  smooth 
except  for  possible  singularities  across  the  characteristic 
cone  in  the  space  z  >  0.  We  expect  that  there  will  be  a 
solution  to  equation  (l6)  which  is  homogeneous  of  degree 
-2  (-5)  for  a+  p+ Y  =  6  (7)  since  5^(x)6^(y  )5'^(z )  is 
homogeneous  of  degree  -(a+p+Y)  and  Q(d/3x,3/Qy,^/3z )  is  a 
homogeneous  operator  of  fourth  degree.   We  therefore  require 
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the  solution  to  have  this  rate  of  decay  at  infinity  except, 
of  course,  along  the  characteristic  cone  in  z  >  0. 

We  will  not  prove  that  the  solution  we  obtain  is 
unique  under  the  above  conditions  although  we  believe  it 
is;  one  difficulty  is  the  dependence  of  the  form  of  the 
equation  on  the  coordinate  system  chosen. 

We  note  that  the  complex  roots  -},      ^(n  )  of  the 
characteristic  polynomial  occur  because  the  equilibrium 
velocity  n^   is  subsonic  with  respect  to  the  fast  distur- 
bance speed  of  the  time  dependent  problem.   We  therefore 
do  not  expect  the  solution  of  equation  (l6)  to  vanish  for 
z  <  0,  as  would  occur  for  equilibria  in  region  I.   We  will, 
however,  not  need  this  as  an  additional  condition  on  our 
solution. 

It  will  be  convenient  for  us  to  solve  equation  (l6) 
with  smaller  exponents  on  the  delta  functions  and  then 
differentiate  appropriately  to  find  the  solution.   We  shall 
use  both  a  +p  +Y  =  3,      ct,p,Y  >_  1,   and  ct  =  p  =  Y  =  1. 

B.   The  Method  of  Fourier  Transforms 

The  use  of  Fourier   transforms  for  the  solution  of 
the  above  problem  will  be  found  wanting  in  several  respects. 
Difficulties  arise  both  in  interpreting  the  result,  which 
is  written  in  terms  of  divergent  integrals,  and  in  showing 
the  solution  is  indeed  correct  i.e.  justifying  such  procedures 
as  the  displacement  of  the  contour,  etc.   The  interpretation 
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of  the  divergent  integrals  is  a  serious  difficulty  in  that 
the  correct  interpretation  should  give  a  surface  of 
singularity,  namely  the  surface  formed  by  the  tangent  half 
lines  to  the  line  joining  the  slow  loci  in  the  construction 
of  the  characteristic  cone,  and  this  surface  does  not  appear 
in  the  usual  treatment  of  the  above  divergent  integrals,  i.e. 
the  finite  part  integral  considered  by  various  authors,  e.g. 
Hadamard  [11].   In  fact,  this  last  method  is  normally 
expected  to  represent  an  ordinary  function,  while  we  expect 
a  solution  which  is  homogeneous  of  degree  -2  or  -5,  so  it 
will  not  bs  regular  at  the  origin  and  must  be  expected  to  be 
a  distribution.   The  method  we  shall  use  later  will  give  a 
solution  in  terms  of  derivatives  of  continuous  functions  so 
it  is  a  well  defined  distribution.   In  addition,  our  method 
will  lead  in  a  natural  way  to  a  proof  that  the  solution  is 
indeed  correct.   However,  Fourier  transforms  have  been  so 
often  useful  in  the  past  that  it  will  be  interesting  to  use 
them  here  for  comparisoxi  with  our  later  techniques. 

We  begin  by  seeking  a  solution  v(x,y,z)  to  equation  (l6) 
with  a+p+Y  =  5,  a,^,^  ^  1.   We  transform  this  equation  in  the 
usual  way,  assuming  that  v(x,y,z)  is  such  that  the  boundary 
terms  at  infinity  all  vanish  (even  though  we  expect  that 
the  derivatives  of  v  will  be  singular  across  the  character- 
istic cone  in  z  >  0).   We  solve  for  v(k^,k  ,cu),  the  Fourier 
transform  of  v(x,y,z),  and  the  proposed  solution  is  then 
given  by  the  inverse  transform,  namely: 
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v(x,y,z)  =  Re  ^^-j  JSJ   dk^  dk  dec  e^^^^^^)  ^^^ - 

\       1      _oo  ^  x^  y^  ^ 

with  k^=  (k  ;,k  )  =  kn^j  n''=  (cos  Q,    sin  0),  and  x  =  (x^y)  . 

X  y 

This  integral  is  singular  because  Q,(k,a))  has  the  real  roots 
m  r^^)   =  -k  A-,  „(n^).   For  each  k  ^  0,  we  displace  the  cd 
contour  from  the  real  axis  to  a  contour  c  parallel  to  the 
real  axis  in  the  lower  half  plane  and  above  the  pole  at 
CD  :=  CO  (k^)  =  -k  A^(n*").   This  is  the  proper  displacement  if 
v(xjy,z)  is  to  satisfy  the  causality  condition,  as   will  be 
shown  later.   The  methods  we  use  to  evaluate  the  resulting 
integral  will  follow  closely  those  used  by  H.  Weitzner  for 
the  time  dependent  problem  [8]. 

The  integral  over  cd  may  be  evaluated  by  residues 
to  give  for  z  >  0, 


v(x,y,z)  =  Re  -±A.J:.L  dk  dk  cos"^  9      sin 

Y_l,_>,   ik(5?n"-A  (n^)z) 
A.   (n  )  e        ^ 

k  Q^(iT;-Aj.(n)) 
and  for  z  <  0, 

v(x,y,z)  =  Re  -^— L_  fP  ^^  ^i^  qos"  "^0  sin^"^( 

(27r)^   J  J    X   y 
^  '        -oo 

Y_l^_>^   ik(5f.rf-A^(n^)z) 


Af "(r)  e 

•       ^.        , 


k  Q^(nr-A^(n'')) 
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where  we  have  used  the  homogeniety  of  Q,   to  give 
co.(k'')  =  -k  A.(n'')  and  Q^(lc^-aj)  =  k\^(n^-fDA). 

Considering  first  z  >  0^  we  change  the  variables 
of  integration  to  polar  coordinates  {K,9) ,    -tt  <_  9  <  t,  and 
then  let  k'  =  -k,   0'  =   6+  it     on  (-Tr,0)  to  find 

Y-1  '^   +°° 
v(x,y,z)  =  Re  --^— ^^—  f  d9  f  dk  sgn  (k) 

(2^)     0    -03 

ik(x«n-A . (n  )z) 
cos^~-^0  sin^"^9  ?fj^n^)  _e 


z  >  0. 

The  integral  over  k  exists  as  a  distribution  and 

may  be  evaluated  using 
oo 


,,      /,  N   ikt    2i 
dk  sgn  (kj  e    =  -r-     , 


-00 

[12],  to  give 


v(x,y,z)  =  Re  ±^^-^)  "    J_ 


Y-1    „    cos   0    sm   9  A '.   (n  ) 


de 


{2vf        j=l,2,4  "q     (r.n'-Aj.(n)z)Q^(n',-Aj(n')) 


z  >  0. 


and,  in  the  same  way. 
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p/  ,^y-l    Z           cos^-^Q  sin^-^e  A^"^(n^) 
v(x,y,z)  =  Re  ^^i^-U^   d0  -::^-3^ -^ -—  , 

z  <  0  . 
We  may  rewrite  the  above  integrals  as  integrals 
over  the  unit  circle  by  using  the  relations  A.  {Q-'^) 


(equations  (lO)  and  (15))  and  x*  n]  fl_-jj-  =  -^'   n]  g  ,    which 
give 


/  -,\7-l  p      cos    9   sm^   0  A '.   (n  ) 

v(x,y,z)  =  Re  1--:LZ__  J  ^0  _  X__   "-^.^ -J^ Z^ ^^ 


27r)^    0   J-1^2,4  (x-n-A^.(n  )z)Q^(n, -A^.  (n  )) 

z  >  0, 
(17) 

,    ^.T     f  C0S-I9   sinP-l0   A^-l(r) 

v(x,y,z)  =  Re  ^^^    d©  ^-.^ ::^ -^—   , 

{2Trf     I  (?.n-A^(n^)z)Q^(K;-A^(n^)) 

z  <  0  . 

The  conjugation  signs  do  not  appear  as  we  take  the  real  part 
We  have  the  proposed  solution  expressed  in  equation 

(17)  in  terms  of  divergent  integrals.   We  shall,  by  using 

a  different  method  of  solution,  find  the  proper  interpreta- 
tion of  these  integrals  to  be  that  given  by  formally 
integrating  with  respect  to  x,  y,  and  z,  as  follows!   We 
introduce  the  notation 

P     (x-n-A . (n  )z )  _ 

(18)  P(m,j)  ^        d0 -J -—  log  GAt'-n   -   A.(n^)z) 

0    Q^(nr-A.(n"))        J        J 

where  m  is  a  positive  integer,  j  =  1,2,5,4,  and 
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> 


X    for  j    =   1^2 
G.(x) 
"3         I  X   for  j  =  5j4 


The  solution  is  then 

■1 

r         ^   . 

(27r)^  j=l,2,4 


^  a-1  ^  p-1  ^  T-1   . 

v(x,y,z)  =  Re  (|-)    (^)    (^)    [-^     ZiPd.j)], 
ox     ay     '-'-^      I'p-rr^'^  i=l  .P.li 


z  >  0 
(19) 


.  a-1  .  6-1  .  Y- 


1 


1 


v(x,y,z)  =  -  Re(^)    (^)    {^)         [-^-^   P(l,3)], 

^"^^  z  <  0. 

We  shall  discuss  this  solution  at  length  later, 
giving  the  branches  for  the  logarithms.   In  particular, 
we  shall  show  that  the  functions  P(l,l) +  P(l,2) ,  P(l,3), 
and  P(l,4)  are  continuous  in  x,  y,   and  z, and  thus  v(x,y,z) 
is  a  well  defined  distribution.   The  solution  u(x,y,z)  of 
equation  (l6)  is  then  found  by  taking  the  appropriate 
derivatives  of  v(x_,y,z). 

We  conclude  this  discussion  of  the  Fourier  transform 
method  by  observing  that  if  we  had  displaced  the  co  contour 
in  the  above  derivation  into  the  upper  half  plane,  the 
proposed  solution  v  would  be 


-u  =  ^^  -b  (|?)    (^)    (^)     ^(1^^)'         ^^° 
and 


5  a-1  >  p-1  ^  T-1 
(2^)"=   5^     ^  Sz'    J=l72,5 


-u  =  -  «- T^Ts  <^'   '4^'   '57)   .„?=/(^'J)'  ^'°- 


The  inclusion  of  the  sum  over  the  real  roots  for 
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z  <  0  will  cause  singularities  to  appear  across  the 
characteristic  cone  in  that  half  space,  as  will  be  shown 
later.   This  violates  the  causality  condition  and  this 
otherwise  legitimate  solution  must  be  rejected. 

C.   The  Method  of  Solution  and  an  Example 

¥e  shall  find  a  solution  to  equation  (l6)  by 
extending  the  plane  wave  representation  for  the  fundamental 
solution  of  homogeneous  hyperbolic  equations,  as  given  by 
F.  John  [15],  to  a  representation  for  the  solution  of 
equation  (I6)  with  ct  =  p  =  7  =  1. 

We  begin  by  assuming  for  the  moment  that  equation  (I6) 
is  hyperbolic  i.e.  that  all  the  roots  A.  (n^)  of  Q,(n^?N)  are 
real.   The  solution  to  equation  (I6)  with  a  =  p  =  Y  =  1, 
using  Duhamel's  principle,  is  H(z  )u-,  (x,y,  z) ,  where  H(z)  is 
the  Heaviside  function  and  u-,(x,y,z)  is  a  solution  of  the 
following  Gauchy  problem  in  the  space  z  >  0! 


'fe  '  S^'   ^)ui(x.y,z)  =  0  ,         z  >  0 


(20)  with  initial  conditions 


a^S^u.  (x^y^O"^)      fo         for   i  =  0,1,2 


Sz^  |5(x)5(y)   for  1    =   ^ 

The  solution  u-,  of  this  problem  is  given  by 

F.  John  [13]  in  terms  of  the  two  dimensional  Laplacian 

2     2 
(A_  =  D  +  D  )  iterated  three  times   acting  upon  a 

X     ^     y 

superposition  of  plane  wave  functions,  i.e.  functions 


of  the  form  g(x,n-A.(n  )z).   In  particular^  for  the  case 
of  two  space  dimensions  (such  as  here  with  z  regarded  as 
the  time  variable )j  the  plane  wave  functions  are 
(x^  n^-A. (n  )z)   log  jx  •  n  -  A.(n  )z]  ,   where  m  is  a 
positive  integer.   For  A. (n  )  complex,  these  functions 
are  no  longer  plane  wave  functions  [since  lx^«n  -A.(n^)zl 
=  (x  .n  -A. (n  )z)(x  .n  -A. (n  )z)]   and  u.  is  no  longer  a 
solution  to  the  homogeneous  equation  in  z  >  0.   We  correct 
this  "by  replacing  log  jx^'n  -A.(n'')z|'  by  the  plane  wave 
functions  log  (x"  -n  -A.  (n  )z)  for  i  =  '^,^,    where  we  use 
the  fact  that  the  roots  -A.  (n  )  of  the  polynomial  Q,(n^,A) 
are  analytic  functions  of  9    to  extend  the  logarithms  into 
the  complex  plane.   We  will,  of  course,  have  to  give  branches 
for  the  logarithms.   We  make  this  change  in  u-,(x,y,z),  take 
the  real  part,  and  arrive  at  a  new  function  Up,  where 
Q,(S/^x,  S/^y,  S/5z)u„  =  0  in  z  >  0  (with  A^  ii(n^)  complex). 

The  function  u  unfortunately  does  not  satisfy  the 
initial  conditions  of  equation  (20).   In  fact,  we  will  have 

aQ^\^{x,y,0^)  rf.(x,y)        for  1=0,1,2 

bz^  |f^(x,y)+5(x)5(y)   for  1=3, 

with  f^  ^  0,   i  =  1,2,3,4.   The  function  H(z)u2  and  its 

first  three  normal  derivatives  will  then  not  be  continuous 

2     2 
across  z  =  0,  x   +  y  7^  0,  and  H(z)up  will  not  satisfy 

equation  (I6).   We  correct  this  by  finding  a  function 

M(x,y,z)  which  is  a  solution  to  the  homogeneous  equation 


■hh- 


for  z  ^  0  and  which  will  smooth  H(z)up,  i.e.  so  that 

u^  =  H(z)u2  +  M  will  have  at  least  three  continuous  normal 

2    2 
derivatives  across  z  =  0^  x  +  y  ^  0.  The   function  M 

may  be  found  by  assuming  it  to  be  a  superposition  of  plane 

wave  functions  similar  in  form  to  Up^  but  using  only  the 

plane  waves  x-  n  -  A-,  ^(njz  in  order  not  to  violate  the 

causality  condition.   The  arbitrary  functions  included  in 

this  superposition  will  be  determined  so  that  Up  and  its 

first  three  normal  derivatives  are  continuous  across  z  =  0, 

2    2 

X   +y  y"  0.   The  function  u^  will  then  have  the  property 


3"    ,     .;,    a^uj 


that  its  jumps,  ^^  u  (x,y;,0  )  -   ^^  (x,y,0  ),   across  z  =  0 

vanish  for  n  =  0,1,2,  and  equal  —  6(x)6(y)  for  n  =  3. 

^0 

We  can  then  easily  prove  that  u^  is  a  solution  to  equation 

(16)  which  satisfies  the  smoothness,  causality,  and  rate 

of  decay  at  infinity  conditions  given  above. 

We  illustrate  this  method  by  using  it  to  find  the 

fundamental  solution  to  Laplace's  equation  in  three 

dimensions,  i.e.  a  function  u(x,y,z)  such  that 

u   +  u   +  u   =  6(x)5(y)5(z) .   The  characteristic 
'^^        yy    2,z 

polynomial  is  Q,(n, A)  =  A~  +  1,  with  roots  -A  (n^)  =  +i. 
If  the  roots  were  real,  the  fundamental  solution  would  be 

H(z)u,   where,  with  A   =  D  +  D   , 

1  2   x^    ^    y 

A      Oil-       ,—>   — >  ^   vj-,    I— >  ^*,    I 
-=*   ^"       (x  •  n  -A.zj-^log  x  •  n-A,z 

+  + 


H(z)u^  =  H(z)  —-^-2  J  d0  [ 


3i(27r)^  "q  2A_^ 


{^'n'-l\   z)^logl  }?"•  rT -A  z] 

+  ~ ~—   ] 

2A 
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With  -A  =  +  ij  H(z)u-,  is  not  a  solution  of  the  homogeneous 

equation  for  z  >  0,  so  we  extend  the  logarithms  into  the 

complex  plane  to  find  the  new  function 
H(z)  A^   27r 
H(z)u^  =   Re ^ —  f  d0  [  (x^n  +  iz)-^log(x".rf+ iz) 

-  (x  •  n  -iz)-^  log  (x«n  -iz)]  . 


The  logarithms  are  defined  to  be  the  principal  branches 
cut  on  the  negative  real  axis  i.e.  -ir   <   Im  log  z  <_  ir.   This 
function  is  not  smooth  across  the  plane  z  ==  0^  x  +  y  ^  0. 
We  then  seek  a  function  M  which  is  a  homogeneous  solution 

for  z  y  0  and  such  that  u-,  ^  H(z)u„  +  M  and  its  first  three 

2    2 
normal  derivatives  are  continuous  across  z  =  0^  x  +  y  =0. 

We  assume  M  has  the  form 
2ir 
M  =  Re  A_^  J  d9  [a(9)(x-n  +iz)-^log(x'n  +iz) 


^   0 


+  b(9)  (5r.rr-iz)^log  (x^'rT-iz)] 


with  the  logarithms  defined  as  above.   This  function  is  a 
solution  of  the  homogeneous  equation  for  z  ^  0  since  it  is 
a  superposition  of  plane  waves.   We  now  pick  a (9)  and  b(9) 
so  that  H(z)u„  +  M  and  its  first  three  normal  derivatives 
are  continuous  across  z  =  0.   A  method  for  finding  such 
functions  will  be  given  when  we  consider  equation  (l6)  and 
we  don't  wish  to  enter  the  details  in  this  example.   We 
will  thus  arrive  at  the  function 
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2 

A      27r 


M  =  -  Re -^ 6.9    [  (3?n^+zi)-^log(x'«rf +zi) 

2(27r)^3:2i  Q 

-  (x'.n'-zi)^  log  (x*"'  rT -zi)  ]  . 

The  proposed  solution  is  u  =  H(z)up  +  M  which  we 
may  evaluate  by  applying  the  Laplacians  under  the  integral 
sign  for  z  /  0.   We  then  use  the  identity, 
2v 


6_0_ _    2Fi  sgn  (b) 


0  ^°"  ^  +  ^i         yT+b2 

which  is  true  for  b  real  and  nonzero,  to  find 


^3 


^^./„2,..2,„2 


X  +y  +z 

This,  as  is  well  known,  is  the  correct  answer. 
In  the  case  of  equation  (16),  we  shall  have  to  prove  that 
the  solution  we  obtain  using  this  method  is  correct. 

D.   Derivation  of  the  Solution 

We  apply  the  above  method  to  equation  (l6)  with 
a  =  p  =  Y  =  1.   To  begin,  the  solution  H(z)u-,  to  the 
hyperbolic  problem  is  modified  by  extending  the  logarithms 
involving  the  complex  roots  A^  j.{d)    into  the  complex  plane. 
We  then  show  that  the  resulting  function  H(z)up  is  well 
defined,  using  essentially  the  same  method  employed  by 
F.  John  [13]  to  show  that  H(z)u-^  is  well  defined.   We  find 
the  "smoothing"  function  M,  and  u^  =  H(z)up  +  M  will  be  the 
proposed  solution.   The  proposed  solution  u^  is  shown  to  be 
the  same  function  as  that  fo^ond  formally  by  the  Fourier 
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transform  method  (equation  (19)).   Finally^  we  check  that 
we  have  not  neglected  another  solution  during  one  of  the 
steps  in  the  derivation  of  M.   (This  is  not  a  uniqueness 
proof. ) 

We  begin  with  the  function  H(z)u-,  which  would  be 
the  solution  if  all  the  roots  A. (n  ),  i  =  1,2^3,4  were 
real.   We  take  this  function,  as  given  by  F.  John  [IJ], 
extend  the  logarithms  log  Ix-n  -  A^  K(n)z|   into  the 
complex  plane,  and  take  the  real  part  of  the  resulting 
expression  to  find 

(21)    H(z)u^(x,y,z)  =  -J^^l-Re  ^{     YZP{7.3) 


^  (2Tr)^7:     X   j=l 

where  P(m,j)  was  introduced  above  (equation  (18)).   Branches 
for  the  logarithms  in  P(7j3)  a^d  FiJ ,h)    are  given  by  setting 
log  1=0,  and  placing  the  branch  cut  for  log(x*n  -A^(n  )z) 
in  the  closed  upper  half  plane  and  for  log  (x  "n  -Ah(n  )z) 
in  the  closed  lower  half  plane.  The  logarithms  will  then  be 
continuous  for  z  >  0  since  Im  A^(n  )  =  -  Im  A2,(n  )  >  0. 
For  def initeness,  we  take  both  branch  cuts  along  the 
negative  real  axis. 

We  next  show  that  Up(x,y,z)  is  well  defined  despite 
the  double  points  of  the  real  roots  at  sin  9  =  0  at  which 
points  the  denominators  of  PiYfl)    and  P(7j2)  vanish.   Each 
of  these  integrals  is  singular^  however,  their  sum 
P(7,l)  +  P(7,2)  or,  more  precisely,  S,  the  sum  of  their 
integrands,  is  continuous  and  we  shall  always  use  these 
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Integrals  in  this  combination. 

The  GOxatinuity  (and  differentiability  for  m  >  2)  of 
the  sum  may  be  shown,  follovjing  F.  John  [15],  more  generally  for 
P(m_,l)  +  P(m,2)  with  m  _>  2.   We  first  note  that,  for 
sin  0  /  0,  a  residue  expansion  gives 


>m 


S  =  1 


(x  •  n-A.(n  )z)  log(x  •  n-A.(n  )z) 


f  — >  ,  /  — ^>  • 


(22)  =    -  ^ 


1=1,2  Q^(n,-A^(n  )) 

(x^*  n-A.  (n^)z  )'^log(x  •  n-A.(n^)z) 


i=^A  Q^(nr-A^(n'')) 

+  _1_  r  dA  (^"^  n"-Az)"^log(5f-  n-Az) 


— >  — > 

X  •  n 


The  path  c  is  a  circle  of  radius  R  about  the  point  7\^= 

slit  along  the  parallel  to  the  positive  imaginary  axis  through 

A^  for  Re  A^  ^  }\^     and  along  some  other  radius,  say 

arg  (A-Ap^)  =  77-  ,  for  Re  A^  =  A^^  ,  or  c  may  be  any  equivalent 

contour.   R  should  be  so  large  that  c  contains  all  roots 

of  Q(n,A),  e.g. 

R  =  |Aq|  +  2  max  |A^(n)|  ,     i  =  1,2,3,^. 

Along  the  slit  (x-n  -  Az)"^  =  (A  -  Aq)'^z"^  while  Q(n,A) 
vanishes  at  most  like  (A  -  A^-.)  ,  when  the  real  roots  coincide 
at  A^  ,  i.e.  for  A^^  =  0.   The  integral  in  equation  (6)  is 
thus  continuous  in  x,y,  and  z,  since  m  ^  2.   The  sum  over 
the  complex  roots  on  the  right  side  of  equation  (22)  is 
infinitely  diff erentiable  in  x,y,  and  z,  since  the 
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denominators  are  "both  nonzero  and  A^  ^(n  )  are  nonreal. 
The  function  S  is  thus  continuous.   In  order  to  form 
derivatives  with  m  >  2^,  we  write  the  integral  in  equation 
(22)  as 

z""     r  ^^^   a3'^log(a^z)  _  ^m   p    __^_a^ 

'^'  Io^iIr   Q(nra>fA,)       i^        Q(nra>fAo) 

for  Re  'K     =/  },      ,    and  Integrate  the  second  integral  above 

over  arg  (A-A^)  =  jr     for  Re  ^     =  Aq.   The  partial 

derivatives  with  respect  to  x,  y,  and  z  of  these  integrals 

— >  — > 
(recalling  that  Ap,  =  ^^-~—  )  may  be  computed  by  differenti- 

ating  under  the  integral  signs.   If  we  then  integrate  by 

partSj  the  resulting  expressions  will  not  contain  any 

contributions  from  the  limits  of  integration. 

We  thus  see  that  the  integral  in  equation  (22),  and 

hence  S  and  P(m,l)  +  V{vx,2) ,    has  in  the  spi^ce  z  >  0  all 

derivatives  with  respect  to  x,  y,  and  z  of  order  <_  m-2 

and  these  derivatives  are  continuous. 

We  can  moreover  show  that  these  derivatives  have 

4-        2    2 
finite  limits  as  z  — >  0  ,  for  x  +  y  -/  0.   This  may  be  seen 

by  noting  that  as  z  — >  0  ,  we  have  ^q  "~^  °°  *   Therefore, 

for  z„  sufficiently  small,  we  may  take  our  path  of 

integration  for  the  integral  in  equation  (22)  to  be  a 

circle  about  all  the  roots  of  Q,  where  this  circle  does 

not  contain  A^^.   The  derivatives  and  the  limit  z  — >  0 

may  then  be  applied  directly  under  the  integral  sign. 
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In  the  remainder  of  this  thesis ^  the  integrand  of 
P(mjl)  +  P(mj2)j  m  _>  2^  is  to  be  regarded  as  defined  as 
in  equation  (22)  . 

The  function   H(z)Up^   equation  (21)^  for  z  >  0 
is  given  in  terms  of  the  Laplacian  three  times  iterated 
applied  to  a  five  time  continuously  diff erentiable  function. 
The  terms  P(7^l)  +  P(7^2)  are  five  times  continuously 
diff erentiable  by  the  argument  above  and  there  is  no 
difficulty  differentiating  P(7,3)  and  P(7,4).   The 
function  H(z)Up  is  then  a  well  defined  distribution  and 
since  it  is  a  superposition  of  plane  wave  functions^  it 
satisfies  the  homogeneous  equation  for  z  >  0.  We  are  using 
the  fact  that  H(z)Up  is  a  distribution  to  interchange  the 
operator  Q  with  the  Laplacians. 

We  next  seek  a  function  M  such  that  M  satisfies 
the  homogeneous  equation  for  z  ^  0  and  such  that  u^  =  H(z)up+M 

and  its  first  three  normal  derivatives  are  continuous  across 

2    2 
z  =  0,  X  +  y  7^  0.   We  assume  M  is  of  the  form 

27r 
M(x,y,z)  = ^^  Re  A   [  r  d0 

•  a(0)(x':  n'--h^{n')z)^   log  (x':  n^-\^{n')z) 
+b(0)  (x^.  rf -Aj^(n^)z) '  log  (x^*  n^-A2j^(n^)z)  ] 

for  z  =/  Oj  where  a  (9)  and  b(9)  are  to  be  determined. 
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¥e      assume,  without  loss  of  generality,  that  a(0)  and 
b(0)  are  continued  in  9   with  period  27r,  i.e.  a(9-f2Tr)  =  a  (9) 
and  b(0+27r)  =\){Q).      We  do  not  include  integrals  over  the 
real  roots,  A-,  ^{n    )    ,    in  M  as  this  will  introduce 
singularities  across  the  characteristic  cone  in  z  <  0, 
violating  the  causality  condition.   We  take  the  same 
branches  for  the  logarithms  here  as  were  taken  for  the 
logarithms  in  \{{z)\x^,    equation  (21).   This  choice  is  made 
to  avoid  introducing  singularities  in  M  for  z  ^  0  due  to 
branch  cuts  and  so  that  all  logarithms  we  consider  will 
be  defined  in  the  same  way. 

The  functions  a((9)  and  b(0)  will  now  be  determined 

by  the  condition  that  H(z)up  +  M  and  its  first  three  z 

2    2 
derivatives  be  continuous  across  z  =  0,  x  +  y  =/  0. 

(We  will  show  later  that  H(z)up  +  M  is  analytic  in  a 

neighborhood  of  z  =  0,  x  +  y  ^  0. )   We  begin  by  finding 

the  limits  of  u^  and  its  first  three  z  derivatives  as 

—   2    2 
z  — >  0  ,  X  +  y   y  0.   The  various  derivatives  and  limits 

of  u^  may  be  taken  under  the  Laplacians  since  u^  is  a 

distribution.   They  may  then  be  taken  under  the  integral 

signs,  since  we  have  shown  the  integrals  appearing  in  Up 

are  five  times  diff erentiable  for  z  >  0  and  these 

+   2    2  / 
derivatives  are  moreover  continuous  asz— >0,x+y  -/O. 

The  integrals  appearing  in  M  may  be  treated  in  the  same  way 

as  they  present  no  difficulties. 

The  limits  as  z  — >  0  of  H(z)u^  and  its  first  three 
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z  derivatives  are  given  by 


n 


4    (-A.( 
log  1 X  •  n  I  2l 


0 

n 


(23) 


i..l  Q^(9,-A^(0)) 


\        tc-^FT^ — vr    d0  (x  .  n  )    (- + 


_>  4.TFEri  J,^^  ^         ^,(^.-^3(9))  Q,(e.-A,(e)) 

Oj<9<_27r 

5     ^n     f     >   >  7  n  ^     (-Ai(e))'' 

+  Re  A  , -5—   f  d0  (x'^  E^)^  "^  ^  ^ ,n=0,1.2,3 

r  {2wfl\      I  i=l   ^^(0,-A.(9)) 

with  Cq  =  0,  c-j_  =  1,  Gg  =  13^  and  c^  =  55. 

Before  we  proceed  in  a  similar  way  with  the 
function  M J  we  will  simplify  the  above  expression.   We  use 
the  following  identity: 


-,q  ^^     m    \}  (   0  for  q  -  0,1, ..m-2 

^  ^   2Tri   J   -QTaT-  -  fzx  ^A^'^i^     I  —    for  q  =  m-1 


where  Q,(a)  is  a  polynomial  of  degree  m  with  leading 
coefficient  a^^  y  0  and  with  all  roots  A.  simple.   The 
path  of  integration  c  is  any  contour  enclosing  all  the 
roots  of  Q.   For  Q  with  coincident  roots,  as  in  our  cage, 
we  define  the  summation  in  equation  (2  4)  by  the  integral 
and  the  identity  still  holds.   Applying  identity  (24)  to 
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equation  (25)  with  n  =  0,1^2,  we  see  that  the  first  and 
third  terms  on  the  right  side  vanish.   With  n  =3,  and  using 
identity  (24),  we  find  these  two  terms  are  the  following! 


^   27r  h, 

->      p     ,— >  — >,4     I— >  — >|   ^3^*  ,—>   — >,4, 
2  "   J  d9  [(x'  n  )   log  X  •  n   +  y,   (x  •  n  )  ] 


aQ4!(2Tr)^  x   ^ 

27r 
(25)  =  ^ ^  A    r^"^  log  r  f  d0  cos^f 


1   A    1  5  X  5  y 

^^^0   X  0 

¥e  note  here  that  if  all  the  roots  A^(0)j  1  =  1,2,3,4, 
were  real,  the  two  terms  we  have  just  considered  would  be 
the  only  terms  on  the  right  in  equation  (23).   The  above 
computation  then  verifies  that  the  function  u-,  given  by 
F.  John  [13]  is  a  solution  of  equation  (20),  and  thus 
H(z)u-,   is  a  solution  of  equation  (I6)   when  the  roots  are 
all  real,  i.e.,  when  the  equation  is  hyperbolic. 

Returning  to  our  case,  we  have  simplified  equation  (23) 
which  now  reads 

S''h(z)u  (x,y,0+)       3    i(-l)^      r     .->  ->^l-rx 

X  •  n  <.0 


(26) 


0<  e<2rr 


.  I 1 + 1 ]  + 

+  -^-   6(x)  5(y)  ,   n  =  0,1,2,3, 
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with  5  ^   =     ]     '         ^  -^ 


i- 


n,3  ■   1.1,  n  =  5 
The  function  M  may  be  treated  in  the  same  way  to  find 

l']M?S^Y_,Jhl  =   ne   ^l ^^^   fdS  [A-5(0)a(9)+A2(9)b(0)] 

Sz^  x''  (2F)^(7-n)l   0 

•  [X '  n    )  log  I  X  •  n  I 


+  Re  A^^  -i-l::!!!^-    J  d9  [-A^(9)a(0)  +  A-2(  e)b  (9 )  ]  (x""-  n^)^"^ 


X  4F(7-n): 
/ov)  X  'n  <0 


0  <_9  <_2rr 


3   =n(-l)"   2'^ 

t 


X   (27r)^7I   ;      ^ 

n  =  0,1,2,3. 


We  now  use  equations  (26)  and  (27)  to  calculate  the 

jumps  of  u^  ■=  H(z)up  +  M  and  its  first  three  z  derivatives 

?    2 
across  z  =  0,  x  +  y  ^  0.   The  equations  that  result  from 

setting  these  jumps  to  zero  are 

Re  A-"   i      d9  X  •  n  )  '  ''  U    \a      ^  /Q\^  -^K-rk — ^  1  a\\ 

X  •  n  <_0  -^ 

0  <_  9  <_2Tr 

28) 


+  2(-A5(9)a(9)  +A2(0)b(9))]  =  0  , 


n  =  0,1,2,3. 
An  immediate  solution  is 
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a(0)  -  -^ 


(29) 

The  proposed  solution  to  equation  (l6)  with  a+  p  +T  =  1  is 
then 

^-> 

u   =  H(z)u   +  M  =   Re ^-  [P(7,  l)+P(7.  2)+  ^(P(7,5)+P(7.4))  ] , 

^        ^  (2Tr)^7i 

(50)  .  ^  ^  0 

A-' 
— > 

u^  =  H(z)u^  +  M  :=  -Re  — Ji^—   [P(7,5)  +  P(7.^)].    z  <  0  . 
-^  2(27r)  71 

This  proposed  solution  appears  different  from  that 

obtained  in  equation  (19).   We  now  show  that  equation  (19) 

and  equation  (50)  represent  the  same  solution  of  equation  (l6). 

A^a-1  rA^^-l  rA^T-l  , 

equation  (19)  occur  because  we  solve  equation  (lb)  with 


First,  the  explicit  derivatives  (;v~)     (^)     (t~)     ^^ 


a-Fp+Y=  5,  a,p,Y>^l   while  here  we  use  a  =  p  =  Y  =  1. 
We  must  therefore  check  equation  (30)  against  equation  (19) 
only  after  dropping  these  derivatives  from  equation  (19) 
(or  adding  them  to  equation  (30)).   Second,  the  functions 
in  brackets  on  the  right  in  equation  (50)  are  five  times 
differentiable  in  x,  y,  and  z,    while  the  corresponding  terms 
in  equation  (19)  are  only  continuous.   This  five-fold 
differentiability  has  been  accomplished  by  using 
rr  ^->   P(^^J)  =    P(I^J)^  J  =  1,2,3,^.   Finally,  we  consider 

X 

the  function  v-,  given  by 
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V 


a' 

— > 

X 


Re  -— -p—  [  YZ     P(7.j)]  ,  z  >  0 


1  "  /^„^2, 


(2Tr)^71   j  =  l,2,4 


and  ^3 


-> 

X 
V 


=  -Re p—  [P(7.3)]  ,  z  <  0  . 

(2-/r)^7: 

The  function  v-,  is  just  the  function  v  modified  by  removing 

the  explicit  derivatives  (d^x)'^"^  (S/^y)^"-^  {^ /dz)^~-^ 

and  by  replacing  P(l,j)  by  ^A   (7,j),  j  =  1.2,3,4. 

'  •   X 
The  difference  between  u^  (equation  (30))and  v-,  is 

— > 

u  -  V  =  Re 1 [P(7,3)  -  P(7,4)].    z  ^  0  , 

which  vanishes  identically  since  P(7,4)  -  P(7,5)  is  pure 
imaginary.   This  may  be  seen  by  using  A|,(9)  =  'h^iOj  ,    and 


Q^Je", -A^(9JT  =  Q^(0.-A|^(9))  .   Our  method  of  solution  there- 
fore gives  the  same  answer  as  that  found  by  a  formal 
integration  of  the  solution  found  by  Fourier  transforms. 

Before  we  proceed  to  check  that  u^  is  a  solution 
to  equation  (l6)  meeting  all  our  conditions,  we  will  show 
that  the  functions  a (9)  and  b(9)  given  in  equation  (29) 
are  the  only  solutions  to  equation  (28),  in  the  sense  that 
any  other  solution  a-,  (0),  b-.(0)  will  lead  to  the  same 
solution  of  equation  (I6).   This  is  not  a  uniqueness  proof 
but  a  check  against  neglecting  some  other  solution  corres- 
ponding to  different  values  of  a(0)  and  b(9). 

We  begin  by  appling  the  Laplacians  in  equation  (28) 
as  follows!   The  integral  in  equation  (28)  is  a  sum  of 
integrals  of  the  form 
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G  =    J  d9  {x^-n^)^'^   g(9)  ,   g(0+2Tr)  =  g(0)  , 

x^  n^<0  n  =  0,1, 2^,5. 

0  <_  9  <_  27r 

We  may  write  this,  using  polar  coordinates,  r  =  |x^U 

,  — > 

if/  =  arg  X,    as 

H  3^/2 
G     =  J     d0    r^~"   cos^"^(e-^)    g(9)    . 

■lf/+TT/2 

The   application   of     A        =    {-^   — ^   +  -  -^  +  — p)        gives 

X  r     ^Tp  Br"^ 

Af_^G     =      (7-n)(7-n-l)(7-n-2)(7-n-5) 

X 

?//+5^/2 
.         J        d9   r''-''"^   cos^"^-^(0-^)    g(9) 
^+7r/2 

The  application  of  a  third  Laplacian  will  give 
end  point  contributions  for  n  =  2,3  and  we  find 

A^^G  --      (6^^Q+  5^^^)  7:    J  d0  rl-'^cosl-^(0-^Og(9) 


(31) 

r^ 

n  .=  0,1,2,3. 


-  5^^2  ^^  [g(^^H|)+g(^+|^)]-  5^^^  it-  [g:(^+l^)-g.(^-  I)], 


We  will  need  to  consider  equation  (2S)  only  for 
n  =  2,3  to  show  that  a (9)  and  b(9)  are  uniquely  determined 
in  the  above  sense.   We  use  equation  (31 )  to  show  equation 
(28),  with  n  =  2,  is  the  difference  equation 
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sin  (0-^)  [ 


^id) 


A4(0) 


v^Ajreyi "  Q^^^d7^\wjj 


+  2[A^(e)a(9)  -A^(9)b(9)]] 


d=Tp+Tr/2 


=      0 


Solving  for  Ar(9)a(0)  -  '}\j^{e)h  {9) ,    we  find  the  general 
solution  of  this  equation  as  a  particular  solution^  given 
by  equation  (29) ^  plus  the  general  solution  of  the 
homogeneous  equation,  which  is  any  function  A (9)  which 
satisfies 


(32) 


A(9+  tt) 


-  A(0) 


The  functions  a (9)  and  b(9)  therefore  satisfy 


2, 


-A^(0)  Ag(9) 


(33)   A^(0)a(9)  -  A^(9)b(0)  =  2Q^^rf7TXJCe-n"  "  2Q^(  9, -A^  (  9  ) )  +^  ^  ^  ^ 

For  n  =  5;,  we  use  equation  (31)  to  show  equation  (28) 
is  the  differential  difference  equation 


d   .     b'(^) 
W    L  Q^(9,-A^(e 


,3 


A^(e) 


■A 


4' 


TT 


+  2(A^(9)a(9)  -  Ai;(9)b(9))] 


0  =  ^+37r/2 
9=^+7r/2 


=  0 


Solving  for  A^(0 )a (9 )  -  A^(9)b (9) ,  we  find  the 
general  solution  as  a  particular  solution,  given  by 
equation  (29),  plus  the  general  solution  of  the  homogeneous 
equation,  which  is  any  function  B(9)  which  satisfies 
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{J>h)  B(0+7r)    ^   B(e)    +  G    ,  c        constant. 

This   gives 

(55)  A5(0)a(9)  -  A5(9)Me)  =  2Q;ji7r,-^j- - -2Q^.-^(-eyy  -©(e) 

This  equation  gives  B(9)  in  terms  of  functions 
all  with  period  2it ,   while  equation  {j)^)    gives 

B(0+2Tr)  =  B(9+7r)  +  c  -  B(e)  +  2c 
so 

c  =  0 
and 

(36)  B(9+  it)  -  B(0)  . 

Equations  (33)  and  (35)  are  linear  equations  for 

a(9)  and  h(9)  with  the  solution 

.  A(e)A.(0)  -B(9) 

a(9)  = i +  -2 

2Q  (9,-A  (0))    A;(9)(A^(e)-A  (9)) 

(37)  A     ^       ;?     H     ;? 

A(e)A  (0)  -B(9) 
2Q^(9.-A^(0))    A^(0)(A^(e)-A^(0)) 

This  solution  is  unique  because  the  determinant  of 

2     2 
the  linear  system  is  -2Ah(9)A^(0)  Im  (A^(9)),   which  is 

never  zero  since  'KhiQ)    =  ^"(S)  is  not  real.  The  particular 

solutiOxn  we  have  selec;ted  in  equation  (29)  corresponds  to 

A(0)  s  B(0)  s  0.   If  we  choose  other  values  of  A(9)  and  B{e) , 

subject  to  equations  (52)  and  {j>6) ,    we  will  find  another 

•X- 

proposed  solution  u  different  from  u  as  given  m 
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equation  (jO),   However,  we  now  show  that  u^  -  u  will  be 
a  polynomial  of  degree  less  than  two.   Hence  both  u  and  u., 
will  lead  to  the  same  solution  u  of  equation  (32) ,  because 
u  is  found  by  applying  three  or  four  (depending  on  a+p+Y) 
derivatives  to  u^ . 

To  see  that  u* -  u  is  a  polynomial  of  degree  less 
than  two,  we  first  collect  those  terms  in  the  difference 
which  contain  A(0),  which  gives 


1 


(u^-u^)  ^  ^  =  — —-f;   Re 
5   5^A(0)    (2Tr)2 


d0  A(9)  2i 


( — >  — >, 


A.(9)(x".n-A^(0)z)log(x.  n-A  (0)z) 


2 


A"(9)(A,(0)-A^(0)) 


3 


^4 


3 


where  we  have  applied  the  Laplacians.   We  rewrite  this 

integral  on  the  interval  (■n',27r)  by  letting  9'  =  0  -'w , 

and  then  use  the  relationships  A(9+Tr)  =  -A{9), 

We    thus    find 


A^^^(0+7r 


■A^^^O) 


and   X  •  n 


9+Tr 


— >  — >  I 

-X  •  n 


T 


(u*-u^) 


Re 


d9   A (9)    2i 


^      ^■a(9)         (27r)' 

A^(9)(3r.n-A^(9)z)    logC?- n'-A^(9  )z  ) 


Im    [Re    [ 


0 


3 


^5(0)(A4(9)-A3(9)) 


->   -> 


A.(9)(^.i!-A    (9)z) 

+  TT    i   -| 2 ] 

A3(9)(A^(9)-A    (9)) 


27r 


(27r) 


^  2 


Re 


d9  A (9)    2iF  Im  i 


A.(9)(5^.f?-A^(9)z) 


3 


^i0){\{9)-}.^{9)) 
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which  is  a  first  degree  polynomial  in  x,  y,  and  z.   The 
same  method  shows  the  terms  involving  B  also  only  contribute 
a  first  degree  polynomial.   We  have  therefore  not  neglected 
any  other  solution  to  equation  (l6)  when  we  chose  a (9)  and 
b(9)  as  in  equation  (29). 

E.   Verification  of  the  Solution. 


We  must  now  show  that  u^  is  a  solution  of  equation  (l6) 


with  a  =  p  =  Y  =  1  and  hence  that 
(38)  u  ^  ^ 


Sx^-l^yP-lBz'^-l    ' 
(where  here  a+p+ Y  =  6  or  7)  is  a  solution  of  equation  (l6) 
and,  further,  that  u  satisfies  the  smoothness,  causality, 
and  decay  at  infinity  conditions  given  in  section  5-A . 

We  show  in  this  subsection  that  u  is  analytic  for  z  <  0 

2    2 
and  in  a  neighborhood  of  z  ^  0,  x  +  y  ^  0,  and  that  u 

satisfies  equation  (l6). 

We  begin  by  showing  that  u  satisfies  the  causality 

condition.   The  function  u^  is  defined  as  three  Laplacians 

acting  upon  a  function  w^  given  by 

w^  -  —V—  Re  (P(7,l)  +P(7,2)+^  (P(7.3)  +P(73^)).   z  >  0 
5   (27r)^71  ^ 

and 
w   =  -  —  -J Re  (P(7,3)  +  P(7,^))  ,  z  <  0. 
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We  first  show  that  w^  has  five  continuous  derivatives 
with  respect  to  x,  y^  and  z  for  z  >  0  and  that  w^  is  analytic 

in  x^  y^  and  z  for  both  z  <  0  and  in  some  neighborhood  of 

2    2    2 
the  plane  z  =  0^  x  +  y  +  z   7^  0.   We  have  already  seen  that 

P(7,l)  +  P(7,2)  has  five  continuous  derivatives  for  z  >  0 

+    22 
which  all  have  limits  as  z  — >  0  ^  x  +  y   ^  0.   The  functions 

P(7^5)  and  V{^ ,k)    are  analytic  in  x^  y^  and  z  for  z  7^  0  since 

their  integrands  are  analytic  in  x,  y  and  z  for  each  0<_9j<27r 

and  these  integrands  are  continuous  in  Q,    x,    y,    and  z,  [1''^]. 

Moreover,  P(7,5)  and  P(7,4)  and  all  their  derivatives  have 

■i   2    2 
limits  as  z  — >  0  ,  x  +  y  ^  0,    and  these  derivatives  and 

limits  may  all  be  taken  under  the  integral  sign. 

2   2 
We  now  consider  the  function  u^  near  z  =  0,  x  +y  /  0. 

We  have  chosen  M  such  that  u^  and  its  first  three  normal 

derivatives  are  continuous  across  z  =  0.   We  will  be  able 


to  apply  the  Cauchy-Kowalewski  theorem  if  we  can  show  that 

S  w^ 
on  z  =  0,  ^; — r^  ,  n  =  0jlj2j5  are  analytic  in  x  and  y. 


2   2  2    2 

x  +y  y  0.   On  z  =  0,  X  +  y  /  0,  the  values  of  these 

functions  may  be  found  by  taking  the  limit  from  below,  i.e. 

as  z  — >  0~,  which  by  equation  (27)  and  identity  (31)  gives, 

in  polar  coordinates. 


S^ 


w 


3 


oz 


Re 


z=0   2(27r)^(7-n)  1 


(39) 


-r 


T-n 


log  r 


f    (-A.(0)) 


n 


0 
^+37r/2 


d0[. 


(-A,(0)) 


n 


+ 


^(0,-A^(9)}-^Q^(0,-A^(0)) 


+  TTir 


7 


-(-^3(0)) 


n 


vH 


d0[o  ,n^     .   /n\V  + 


^+7r/2 


_^_^_^    (-A4(0)) 

n=0,l,2,3. 
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These  functions  are  all  analytic  in  r  and  f   for  r  ^0, 
since  the  indefinite  integral  of  an  analytic  function^  as 
above,  is  analytic.   The  functiorB  A-,  h(9)  are  analytic 
since  they  are  non-coincident  roots  of  a  polynomial  with 
coefficients  analytic  in  9   and  with  highest  coefficient 
nonzero.   The  plane  z  =  0  is  not  characteristic  for  our 
equation  (this  was  a  condition  in  choosing  our  coordinate 
system,  and  in  particular,  this  follows  from  the  condition 
that  a^   never  vanish,  so  all  roots  A. (9)  are  finite). 
¥e  may  therefore  apply  the  Gauchy-KowalQAjski  theorem,  which 

says  that  there  is  a  unique  analytic  solution  of  equation 

2    2 
(16)  with  a+p+Y  =  1  in  the  neighborhood  ofz  =  0,  x+y  =/  0 

2    2 
with  data  onz=0,  x+y  /O  given  by  three  Laplacians 

applied  to  the  functions  given  in  equation  (39).   This  then 
must  be  our  solution  u^ . 

The  above  arguments  show  that  u^  and  hence  u  is 
analytic  for  z  <  0  and  thus  u  satisfies  the  causality  condi- 
tion.  We  have  actually  shown  more,  i.e.  that  u  is  also 

2.   2 


ana 


lytic  in  a  neighborhood  ofz  =  0,  x+y  =^  0   and  that 


u^  is  five  times  continuously  diff erentiable  for  z  >  0. 
¥e  will  now  use  this  information  to  show  that  u^  is  a 
solution  of  equation  (l6)  with  a  =  p  =  Y  =  1  and  hence 
u  is  a  solution  with  a+p>+Y  =6  or  7,  a,p,Y  ^  1. 

We  will  need  some  preliminary  results  on  integrals 
of  the  form 
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+00  +00  >N 

(^0)1=  J   J  dx  dy  (A_g(x,y,0))  [       2.    ] 

-GO  -co  ^  °^  °^  ^^ 

where  0  _<  N  =  a+b+c  <_  '^ ,      g(xjy,z)  is  any  C?^    function 
(i.e.  infinitely  diff erentiable  with  compact  support)  and 
[fJ^^Q  =  f(xjy,0  )  -  f{x,y,0~)      is  the  jump  in  f  across 
z  =  0.   We  claim  all  integrals  of  the  form  I  are  zero 
except  when  the  jump  function  is  [^'^w^/dz^]  _„  where 


the  integral  equals  —  g(O^OjO).  We  may  calculate  the 
jump  functions  in  the  same  manner  we  used  for  equations 
(26)  and  (27),  to  find 

(^1)  [-4~b-ir]       =  «-=  --2  [  4^TT^-^  +  S-^""! 


0 


de  cos-e  sin"e  g  ^^(e:.,.(9)) 


which  is  identically  zero  by  equation  (24)  for  c  =  0,1,2, 
The  only  nonvanishing  jump  function  is  then 

b\  r'^  lo  r      4 

^^h=0   =  f^lf^  +  D  r   ,     D  constant 

by  equations  (4l)  and  (24). 

We  may  now  evaluate  the  only  nonzero  integral  of 
form  I,  which  we  will  call  K,  to  be 
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X  bz^    'z=0 


K  =  J     J  dx   dy    (A^^g(x,y,0))     [-^] 


(^2) 


-co  -oo 


+00      +00 

=  I     I  dx   dy   A_^g(x,y,0)[2:|^log   r   +  D^  ] 
-oo    -00 


;(o.o,o) 


Dt       constant. 


¥e  have  performed  several  integrations  by  parts  and  used 
that  (log  r)/27r  is  the  fundamental  solution  of  Laplace's 
equation   in  two  dimensions. 

The  proof  that  u^  is  a  solution  of  equation  (52) 
with  a  =  [3  =  Y=lis  now  straightforward.  We  consider 
the  following   integral  M  over  the   entire    (x^y^z)    space  V, 


-  SH  f(^.y^^)  Q(|^  ^  I-  >  h  ^3^-^y^-)  dv 


M  =  M  I  f(x,y,z)  Q(^  ,  y-    ^ 
V 


oo 
■where  f  is  any  C„   function.   This,  hy  the  definition  of 

the  derivatives  of  a  distribution ,  is 

^-Ul    [^->Q(|r^^^^)  f(x,y,z)]  w^(x,y,z)  dV . 
V     "" 

We  now  write  M  as  M  =  M^  +  M^  where  M  is  the 

integral  over  z  >  0  and  M  is  over  z  <  0.   We  have  shown 

that  w^  is  at  least  five  times  continuously  diff erentiable 

for  z  -/  Oj  so  we  may  integrate  by  parts  repeatedly  in  M 

and  M  to  cast  the  operator  Q(^/dz,  S/dy,  d/5z)  on  the 

function  w^.   We  will  be  left  first  with  the  volume 

integrals 
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N     ^  -     ---'         '         ' 


=  111   [A,  f(x.y,z)]  Q(|^.^.^)W3  dv=  0  , 


z  >  0 


z  <  0 
since  w^  is  a  superposition  of  plane  waves  and,  therefore, 

a  solution  to  the  homogeneous  equation  for  z  7^  0.   We  will 

also  have  many  boundary  terms.   Those  from  infinity  vanish 

because  f  has  compact  support.   The  remaining  boundary  terms 

can  be  grouped  in  pairs,  one  mexmber  of  each  pair  arising 

from  M  and  the  other  from  M  ,  and  each  pair  may  be  written 

in  the  form  I,  equation  (4o),  with  g  as  some  constant  times 

some  derivative  of  f.   These  will  then  all  be  zero,  by  our 

argument  above,  except  for  those  pairs  which  are  of  the 

form  K.   Such  a  pair  can,  however,  arise  only  from  the 

term  in  Q,(S/3x,  S/^y,  d/5z)  proportional  to  B  /dz    .      VJe 

then  have  as  the  only  nonzero  boundary  term 

aQ  JJ  dx  dy  (A£^f(x,y,0))[-^]     =  f (0,0,0) 


X 


"Sz-^  z=0 


by  equation  (^2) . 

We  thus  have  shown  that 


M  =  JJJ   f  (x,y,z)  Q(3^,  5^  ^57)  ^^(^.y.z)  dV  =  f  (0,0,0) 
V 

On 
for  any  Gq  function  f.   This,  however,  is  the  definition 

of  equation  ^l6)  with  a  =  p  =  Y  =  1,  for  u^  a  distribution. 

We  thus  have  shown  that  u^  is  a  solutioxa  of  equation  (I6) 

with  a  =  p  =  Y  =  1  and  that  u  =  S^^+^+^-^u^/  b^-\^^-\  ^^'^z 

is   a    solution   of    (I6)   with  a+p+Y=6    or  7,  a,P,Y  ^  !• 
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F.   The  Singularity  Surfaces. 

The  function  u  has  now  been  shown  to  satisfy  equation 
(16)  and  the  causality  condition.   In  addition,  u  is 
homogeneous  of  degree  -2  (-5)  for  a+  p +  Y  =  6  (7)  so  it 
will  have  the  specified  rate  of  decay  at  infinity  along  any 
path  on  which  it  is  continuous.   We  are  then  left  with 
showing  that  u  is  smooth  except  on  the  characteristic  cone 
in  z  >  0.   This  will  complete  our  proof  that  u  is  a  solution 
to  equation  (I6)  which  satisfies  all  the  conditions  set  forth 
in  Section  5-A. 

We  start  "by  showing  that  u  is  smooth  in  all  of  space 
except  possibly  for  those  points  which  do  not  meet  both  of 
two  conditions.   The  locus  of  points  not  satisfying  both 
conditions  is  then  given  by  geometric  constructions.   We 
next  regard  the  characteristic  cone,  given  by  our  construction 
(Section  3),  in  the  space  z  >  0  as  the  sum  of  two  surfaces, 
S-,  and  S„,  and  we  show  that  those  points  that  do  not  satisfy 
our  first  condition  from  the  surface  S,  and  those  not 
satisfying  the  second  condition  form  the  surface  Sp.   The 
behavior  of  u  near  these  surfaces  will  be  discussed  in  the 
next  section. 

We  have  already  seen  that  u  is  analytic  for  z  <  0 

2    2  • 
and  in  a  neighborhood  of  the  plane  z  =  0,  x  +  y  "/  0. 

Moreover,  we  also  know  that  the  functions  P(7:>3)  and  P(7;>'^) 

are  analytic  in  x,  y  and  z  for  z  ^  0.   Therefore,  the  only 

terms  in  u  which  may  be  singular  are,  by  equations  (58)  and 
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(30),  the  terms, 

(^3)   J  s  Re  ^^-^        -^^     ^  ^  [P(7,l)+  P(7,2)] 

Sx"    dyP  ^  ^z"^     X   (27r)^7: 

with  z  >  0.   (We  are  excluding  in  this  discussion  the  origin 
and  the  singularity  concentrated  there.) 

Let  us  begin  the  description  of  the  singularity 
surfaces  of  J  "by  looking  more  closely  at  the  function 

^t^  {x"'   E-A.  (0)z)''loglx':  F-A.  (9)z 

(^4)   P(7.l)+  P(7.2)  =  r  d0  H-        ^  ^ 


0 


i=1.2        Q^(9,-A^(9)) 


We  consider  a  fixed  point  (x^,y^,Zp^)  and  we  attempt 
to  discover  whether  P(7jl)  +  P(7,2)  may  be  differentiated 
infinitely  often  in  some  sufficiently  small  neighborhood  N 
of  x„.   Differentiating  repeatedly  under  the  integral  sign 
in  equation  (4^)  will  not  suffice  to  prove  this  since  we 
have  already  seen  that  we  can  guarantee  only  five  continuous 
derivatives  of  the  integrand.   The  difficulty  is  the  presence 
of  the  double  points  at  sin  0=0,  where  A-,  p  =  0.   Repeated 
differentiations  under  the  integral  sign  will  also  cause 
singularities  to  appear  at  the  zeros,  if  any,  of  the  f-unctions 
f^(x^0)  5  x":   n-A^(9)z,  i  =  1,2,  0   <   9    <   2Tr .      We  call  the  roots 
of  f .  (x^0)  for  0  _<  0  <  27r  singular  points  and  denote  them 
by  0^.,  i  =  1,2,  and  j  =  l,2,...,m^. 

We  will  be  able  to  differentiate  P(7,l)  +  P(7,2) 
repeatedly  after  moving  the  contour  of  integration  off  the 
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real  axis.   To  do  this^,  we  extend  the  logarithms  into  the 

— >  — > 


complex  plane,  i.e.  we  replace  log  [x  •n-A.(0)z|  by 
log  (x  •  n-A^(0)z)  with  suitably  defined  branches.   This 
procedure  will  succeed  if  two  conditions  are  met  and  the 
points  which  do  not  satisfy  these  conditions  will  form  the 
surfaces  of  singularity  for  P(7^l)  +  P(7^2)  and  hence  for  u. 

The  first  condition  arises  when  we  define  the  functions 
log  If^ix'^e)],    i  =  1,2,  for  9   complex  and  for  all  points  x^ 
in  N.   We  define  the  logarithms  by  taking 
-  7r/2  <  arg  log  f^(xQ,0)  <   J>ir/2,    i  =  1,2,  and  the 
logarithms  for  all  other  x  e  N  are  then  given  by  analytic 
continuation.   This  definition  assures  that  the  logarithms 
coincide  at  the  double  points  Q   =    0,it,    for  x  =  x^^,  as  they 
must  if  the  integrand  of  P(7,l)  +  P(7^2)  is  not  to  be 
singular.   Our  first  condition  is  that  the  logaritlims 
coincide  for  all  points  x  e  N  when  0  =  0  or  tt. 

Before  we  move  the  9-contour,  we  give  some  results 
about  the  functions   A.  p(9)  and  f.  o^^)*   '^^^  functions 
A-,  q(^)  3^^  analytic  in  G   in  some  small  neighborhood  of 

1,2 

the  real  axis  since  they  are  roots  of  a  polynomial  with 
coefficients  entire  in  9   with  leading  coefficient  never 
zero.   They  are  even  analytic  at  the  double  points  0  =  0,Tr, 
as  we  have  shown  in  Section  4-B.   The  functions f    (x,0) 
then  are  analytic  in  9   in  the  same  neighborhood  of  the 
real  axis.   In  particular,  the  number  of  real  zeros  of  f-,  p 
is  finite  in  any  bounded  region  in  which  f.  p  are  analytic 
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and  thus  the  number  of  singular  points  is  finite.   We  also 
have  that  if  a  is  a  root  of  f^(x^e):,  i  =  1  or  2,  which  is 
sufficiently  close  to  the  real  axis,  then  a  is  also  a  root. 
To  see  this,  we  note  that  the  power  series  expansions  of 
f-,  „(x^0)  about  any  real  9   have  only  real  coefficients, 
since  A-,  ^{9)    and  hence  f-,  ^(0)  are  real  for  real  9. 
This  gives  0  =  f^(x^a)  =  f^(x^a),  i  =  1,2. 

We  now  move  the  9-contour  into  the  complex  plane. 
We  first  translate  the  interval  from  0  ^  9  <_  27r  to 
-a  _<  0  <_  27r-  a,  a  ^  0,    using  the  2Tr  periodicity  of  the 
integrand.   The  constant  a  is  chosen  so  that  -a  and  27r-a 
are  not  singular  points.   They  will  then  not  belong  to  9. . 
for  all  5r  e  N  since  the  0.  .  are  continuous  in  x  . 

This  translation  places  all  the  singular  points  and 
double  points  in  the  interior  of  the  interval  of  integration. 

We  now  move  the  9-contour  to  a  new  contour  c  which 
is  defined  as  follows:   c  runs  from  -a  to  2Tr-a  on  the  real 
axis  except  for  semicircles  of  radius  e  in  the  upper  half 
plane  about  each  of  the  singular  points  and  about  9^0 
and  IT.      We  choose  e  sufficiently  small  so  that  the  contour 
remains  in  the  region  of  analyticity  of  the  integrand,  the 
semicircles  do  not  overlap,  except  obviously  if  9  .  =  0  or  tt 

J,  J 

for  some  i  and  j,  and  e  is  taken  to  be  smaller  than  the 
minimum  distance  from  the  singular  points  and  from  9=0 
to  the  boundary  points  9  =  -a  and  9  =  2Tr-a .   We  may  make 
this  change  of  contour  without  changing  the  value  of 


-71- 


P{Y ,1)  +  P(7,2)  because  the  Integrand  is  continuous.  There 
will  then  be  no  contribution  from  the  semicircles  as  we  let 
their  radius  e  tend  to  zero. 

We  now  impose  our  second  condition  which  is  that 
the  singular  points  9.   .  remain  real  for  all  x  e  N.   This 
will  be  true  if  the  d.   .  are  first  order  zeros  of  f.(x^,0), 
since,  as  we  have  just  seen,  complex  conjugate  roots  must 
merge  to  form  a  real  root. 

For  e  sufficiently  small  and  x  e  N,  we  will  have 
no  singular  points  of  the  integrand  of  P(7,l)  +  P(7,2)  on 
our  new  contour  if  the  above  two  conditions  are  satis  fied. 
We  may  therefore  differentiate  P(7,l)  +  P(7,2)  infinitely 
often  under  the  integral  sign,.   The  function  J,  and  hence  u, 
is  then  infinitely  diff erentiable  at  the  point  Xq. 

We  now  give  a  geometric  description  of  the  locus  of 
points  which  do  not  satisfy  both  of  the  above  conditions. 
The  surfaces  formed  by  this  locus  are  then  the  only  surfaces 
on  which  u  may  be  singular.   It  will  be  convenient  to  work 
in  p"^=  (x/z,  y/z)  space  since  the  functions  f-,  ^{x,9)    and 
the  solution  are  homogeneous  in  x,  y,  and  z.   We  set 

singularity  surfaces  only  occur  for  z  >  0  (neglecting, 

as  always  here,  the  origin),  so  f-,  „  and  g.  ^  have  the  same 

zeros,  9 .  . . 

We  begin  by  finding  the  zeros  of  g.(p",  9).   We  let 
N-,  p  be  the  normal  speed  loci  defined  in  polar  coordinates 
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by  the  equatiOxis  r-,  p  =  A-,  ^i^) ,    where  r  =  jp  [  and 

9   =   arg  p"  .   We  let   p  '  he  the  inverse  of  p"  in  the  unit 

circle  and  let  I,  „  be  the  inverse  normal  speed  loci 

given  by  the  equations  r-,  p  ==  lA-]  p(^)«   The  zeros  of 

S-i  p(p^^)j  fo^  each  fixed  p  ,    are  found  by  considering 

the  line  £       which  passes  through  n^   with  normal  in  the 

direction  of  p' .   If  i   intersects  I,  „  at  the  point 

P 
^  "  ^0^   ^12"  -'-'^1  2^^0^^  ^"^^^   ^0  ^^  ^  ^^°^    °^  ^1  2^P^^^ 

and  all  real  roots  of  g-,  ^{p,9)    are  given  in  this  way. 
This  is  easily  seen  since  the  line  i   has  equation 


2— >  — >v    ^   ,     — >   .   ^,   P 


p  • (r  s  -  p  )    -    0,    where  s   is  the  running  coordinate^ 
and  if  s  :=  ^; — 7-5 — r-  n  {9    ) ,    i  =  1^2^  lies  on  i      ,    then 

p^'n(9Q)  -^j^(9q)  =  g-j_(pr^o)  "  ^'      ^^-^^  ^-'-^^  shows  that 

if  0^  is  a  real  root  of  g,  ^Cp^e),  then  .r ^5 — r  ^^i^r.) 

u  l,2^r-  ^lj2\^0''      ° 

lies  on  £ 

P^ 
The  second  condition  above  is  then  satisfied  if 

the  singular  points  9.  .    are  simple  zeros  of  g,  nCPn)^)  > 

1 J  1 J  ^   u 

which  is  equivalent  to  the  line  £,_     having  only  simple  inter- 

Pb      _^ 
sections  with  I-,  o-   The  locus  of  points  p   such  that  £ 

P^ 

is  tangent  to  I.  form  a  surface  which  we  call  the  i"^^ 

branch  of  the  characteristic  surface  and  denote  by  C.  , 
i  =  1:,2.   These  two  surfaces  may  also  be  found  as  the 
envelopes  of  the  lines  m.  ^  with  equation  n^ {9)  •  {'s-'K.  {9)nJ^'f) 
=   o,    X  =   \,2.,    with  s   as  the  running  coordinate.   The  lines 
m.  Q   intersect  the  normal  speed  locus  N.   making  a  right 
angle  at  the  point  of  intersection  with  n  (9).   This  is 
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the  envelope  of  the  surfaces  g. (p,9)  =  Oj  i  =  1^2. 

If  a  point  "p^  does  not  lie  on  one  of  the  surfaces  C-,  „j 

the  number  of  real  zeros  of  g,  oCp"-'^)  -"-^  given  by  the 

number  of  tangents  which  may  be  drawn  from  p  to  G-,  p. 

If  we  consider  a  segment  of  C-,  p  which  is  smooth  and 

has  nonzero  curvature^  points  sufficiently  close  to  one 

side  (which  we  call  the  plus  side)  will  have  two  more 

real  tangents  than  points  sufficiently  close  to  the  other 

side  (called  the  minus  side).   Points  on  the  minus  side 

will  then  have  two  complex  conjugate  6    .   which  merge  into 

a  second  order  singular  point  for  points  on  C.  and  then 

give  two  distinct  real  9.  .   for  points  on  the  plus  side. 

The  points  "p^  for  which  the  line  £_     passes  through 

a  point  of  intersection  of  I,  and  Ip  form  a  circle.   In 

our  case,  where  the  only  intersection  is  at  infinity 

(when  N-,  and  Np,  which  are  inverse  to  I-^  and  Ip  respectively, 

intersect  at  9  =  0,7r),  this  circle  and  its  inverse  degenerate 

to  the  line  x  =  0.   The  line  x  =  0  intersects  the  curves 

^     ,^Gi   yCi,    ,    dAj_(0) 

C-^   and  Cp  at  the  points   p^_  =  [- ^  ~-J  =  (0,  —^g )    , 

i      G^     C^ 

since  gj_(p^  ,0)  =  -^^   (p^  ,0)  =  0.   We  will  show  in  the 

next  section  that  the  line  segment  in  p   space  given  by  x  =  0, 
dA.(O)  dA.(O) 

min  — Itt, —  <  y  ^  max  — ^^^ —  is  the  singularity  surface 
i=l,2^^^^-   -1=1,2    ^^ 
corresponding  tu  Lnose  points  which  do  not  satisfy  the 

first  condition.   This,  moreover^  is  that  part  of  the 

characteristic  cone  given, by  our  construction  (Section  3)  , 
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which  is  ruled  by  the  tangent  half  lines  to  the  line 
joining  the  slow  loci  of  the  time  dependent  characteristic 
surface.   We  call  this  line  segment  S-,  (it  is  a  plane 
segment  in  {x,y,z)    space)    and  we   have  that  the  points 
which  don't  satisfy  the  first  condition  lie  on  S-,  . 

We  call  the  remainder  of  the  characteristic  cone 
S  .  It  will  consist  just  of  the  surfaces  C-,  and  Cpj  since 
it  has  been  shown  that  surfaces  given  by  tangent  half  lines 
to  the  characteristic  surface  of  the  time  dependent  problem 
(excluding  the  plane  segments  caused  by  the  presence  of 
double  points)  coincide  with  the  surfaces  given  by  the 
envelope  of  the  functions  f{x,Q)    with  respect  to  0   [  8  ]. 

We  have  now  completed  the  proof  that  the  function  u 
is  a  solution  of  equation  (l6)  which  meets  all  the  conditions 
set  forth  in  Section  5-A  and  u  may  now  be  employed  to  find 
a  solution  to  the  steady  flow  problem  past  a  source^  as 
described  earlier. 
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Section  6.   The  Solution  Near  the  Characteristic  Cone. 
A.   The  Solution  Near  the  Planar  Segment 

To  investigate  the  solution  near  the  characteristic 
cone,  we  first  study  those  points  which  do  not  satisfy  the 
first  condition  given  in  Section  5-F.   This  condition  states 
that  at  the  double  points ^  9   =   0,Tr,    the  two  functions 
log  f.{x,9),    1   =  l,2j  should  coincide  for  all  x   in  a 
sufficiently  small  neighborhood^,  Nj,  of  the  point  under 
consideration,  x^^.   Branches  for  the  logarithms  were  given 
at  X  =  x„  and  the  logarithms  were  then  defined  by  analytic 
continuation  at  all  other  x  e  N.   We  show  that  the  points 
not  satisfying  this  condition  form  the  planar  segment. 

It  will  be  helpful  to  rewrite  P(7.l)  +  P(7:.2)  in 
terms  of  integrals  over  0  <_  0  <_  tt  rather  than  0^9^  2T. 
We  have 

2Tr   ^;r.;?_^  fa\.^lir..-r\^.;^. 


(45)  P(7,l)+P(7,2)  =  ^    d0 


x-n-A^(9)z)'loglx-n-A^(9)z! 


1=1,2  I  Q^(9.-A^(9)) 

^(x.ft-A.(9)z)^logl^.n-A.(e)z| 

0 


1^2%         Q^r--Ajwi 


since  A^(0+7r)  =  -A^(9)  and  Q^(0. -A^  (0+7r) )  =  -q^{9  ,-},^{9)) , 
by  equations  (lO)  and  (I5)  respectively.   As  in  Section  ^-F , 
we  extend  the  logarithms  and  move  the  contour  to  c' 
which  is  that  part  of  c  with  -a  _<  0  <_  ir-a  .   The  only  double 
point  is  now  at  0  =  0  where  A-,  p  vanish. 
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In  this  section^  we  will  always  have  i  =  1^2. 
We  expand  f .  (x^0)  =  x ^«  n ^  -  A^(9)z   in  a  Taylor  series 
about  0  =  0  and.  x^=x^.   Ignoring  second  order  terms^  we 
find 

(46)   f.(xr0)  =  Xq  +  ^  (r-n^-  A.(0)z)l    9  +  (x-Xq)  +  ... 

or 

dA.(O) 
f.  (x,0)  =  X  +  (yQ  +  —^Q—   Zq)©  +  ... 

For  X  sufficiently  close  to  x„  ,/  0  and  9   sufficiently 
small,  we  have  sgn  Re  f .  =  sgn  x  =  sgn  x^.   Therefore  the 
functions  f._,  which  are  the  arguments  of  the  logarithms, 
do  not  cross  the  imaginary  axis  for  x  e  N  and  9   sufficiently 
small.   The  logarithms  will  coincide  at  9  =  0  for  all  x  e  N 
if  they  do  so  at  x  =  x„,  and  the  first  condition  is  satisfied. 

At  x„  =  0,  the  first  condition  will  still  be  satisfied 
if  the  logarithms  still  coincide  at  ©  =  0  after  analytic 
continuation  across  the  plane  x  =  0.   Near  9  -   0,    c'  consists 
of  a  semicircle  of  radius  e  in  the  upper  half  plane  so  that 
Im  0  _>  0.   We  consider  f .  for  x  =  x-,  e  N  and  let  x-,  have 
positive  x-coordinate.   Equation  (46)  shows  that  Re  f.(5c'-,,9)  >0 
and  sgn  (Im  f^(xj,e))  =  sgn(yQ  -  (dA^(0)/d9)  Zq  ).   (We  do  not 


consider  here  the  intersections  of  the  planar  segment  with  G. , 
i.e.  the  points  Xq  =  0,  Yq/zo  =  dA.(0)/d9). 

If  we  now  analytically  continue  log  f.  across  the 
plane  x  =  0  to  a  point  Xp  ,   equation  (46)  shows  the  real  part 
of  f .  is  negative  and  we  again  have 
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dA^(O) 
sgn  (Im  f^{x^,9))    =   sgn  (yq ^ —  z^).   We  no'Aj  let  9   tend 

_>  dAi(O) 

to  zero  and  find  Im  log  t-ix^fO)   =  sgn(yQ-  — g^ —  Zq)^.  Let 

dAi(O)  .    dAi(O) 

yO,max  =  "^f  (-^F-  ^0^  ^^^  ^O^min  =  "^^^  ^^1^"  ^O^- 
If  y  <  y„   .   or  V  >  y       the  logarithms  will  still 

coincide  at  9  =  0  and  the  first  condition  will  be  satisfied, 

On  the  other  hand,  for  y^^^.^  <  y  <  y^^^^^  ,  let 

(P(7^l)  +  P(7,2) )   be  the  analytic  continuation  of 

P(7j1) +  P(7j 2)  across  the  plane  x  =  0.   We  calculate  the 

jump  across  this  plane  obtaining 


(P(7.l)  +  P(7,2))  -  (P(7,l)  +  P(7.2)) 


A 


_>  _>        7         dA.(O) 
-  Re  27riH(-x)  X_    d0 ^ ^ ^ ^ 


J=1.2  ;,  Q^(0.-Aj(0)) 

The  above  integral  is  evaluated  by  moving  c'  onto  the 
real  axis  yielding  -iri  times  the  residue  at  0  =  0  and  a 
principal  value  integral,  which  does  not  contribute  since 
it  is  real.   We  then  have 


:p(7.1)  +  P(7.2))  -  (P(7,l)  +  P(7.2)) 


A 


2^2  ,       .,       .7  dA^(O)  d?v,(0) 

— A-TO) kdoT ;   l^SniYQ-   —^~   ZQ)+sgn(yQ-  ^j^-Zq). 

-hir^   H(-x)(-x)'^ 

dA  (0)     dATToy""        ~: 

The  jump  across  x  =  0  in  the  solution  to  equation  (l6) 


is  then 
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^a+p+T-3  -  H(-x)(-x) 


d9    -    dS-l  1^^,(0) 


^ 


for  y„   .   *=  y^N  <  y^     and  z^  >  0. 

The  planar  segment  is  then  given  by  x  =  0,  y^   .  < v <  y^ 

and  z  >  0. 

dA.(O) 
We  will  now  calculate  the  values  of  — -=-3- —  in  the 

two  coordinate  systems  given  earlier  (Section  4). 

For  the  system  with  z-axis  perpendicular  to  the  direction 

of  A^,    the  values  of  dA.(O)/d0   are  calculated  from 

Q(cos  G,    sin  9,    A)  (equation  (ll).   We  note  that  Q  is 

Independent  of  cos  0. 

Let  u  =  sin  9,    then  in  (u^A)-space  the  line  given  by 

A  =  ptj  u  =  qt,  -00  <  t  <  CO  is  tangent  to   A  =  -A.  (u) 

at  the  origin  if  the  direction  numbers  (p^q)  of  the  line 

satisfy 

(0,0)p2  +  2  Q^,  fO,0)pq  +  Q,,,  fO,0)q^  =  0  , 


where 


,^^^^,^j^       ,  ^  ^^^^^,^^^H  .  ^^^ 


)^^(0,0)  =  -  u^(^  +  -%)  2  sin^*  , 


m    M 

Q   (0,0)  =   Uq(-^  +  ^)  2  sin  *  cos  (t)  , 

m    M 

Q^^(0,0)  =  -2u^(i2  +  ^)    cos^o  +   2 

m    M  m  M 

dA^    dA^    ^ 
^^^^'^  dT-  =  de-  ^U^-9   '   ^^  '^^^^ 
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dA, 


du 


dA. 

1 

u=0 


-  cot  0  + 


?=0 


/(m^+M^)sin^ct) 


yielding  two  distinct  tangents. 

In  the  second  coordinate  system;,  the  same  method 
gives 

P   P  /  '  "  p    p    p   Q        p 

m  M  cos  4)  sin  *  +  /  (m  +M  )m  M  sin  * 

\~   )  ,    —  P   p     p  p    p  • 

^  -  (m  +M  )  -  m  M  cos  4) 

This  coordinate  system  is  only  used,  as  noted  in 
Section  4-B,  when  l-(m  +M  )  +  m  M  cos  *  <  0. 
This  implies  that  the  denominator  above  is  greater   than 
one,  hence  again  giving  two  distinct  tangents. 

B.   Steady  Flow  Across  the  Planar  Segment. 

The  solutions  of  equation  (6)  are  found  by  forming 
the  appropriate  linear  combinations  of  u  for  various  values 
of  a,    p,  Y,  with  a+p+Y  =  6  or  7.   Equation  (47)  shows 
that  unless  ^  =   -y  =   1,    the  jump  u-u.  across  the  planar 
segment  is  zero.  vie   now  show  that  p  >  1  and /or  Y  ^  1 
for  every  term  of  these  linear  combinations,  i.e.  each  term 
is  differentiated  at  least  once  with  respect  to  either  y  or  z 
Thus  the  steady  flow  solution  with  all  its  derivatives  is 
continuous  across  the  planar  segment  of  the  characteristic 
cone. 

For  the  second  coordinate  system, 
Aq  -  AQ(0,sin  <i>,    cos  <i>)    =   AqIIq  ,  and  u^  =  Uq(0,0,1). 
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2     2     2 
The  equation  for  v-,  =  p  (with  A  =  D^  +  D,  ^  +  D^ )  is 

±  A     y     z 


=  P'  -  ^2  AjD^p^  +  ^  AD^(K^.  B^^)  +  -^  A(nJ.V)(nJ.  M^ 


Every  term  on  the  right  side  is  differentiated  at 
least  once  with  respect  to  either  y  or  z  giving  the  result 
for  V-,  =  p. 

Equation  (5A)  gives  v„  =  a  in  terms  of  p   and  p  , 
so  a  is  also  smooth  across  the  planar  segment.   The  equation 
satisfied  by  v^  =  p  is 


[D^  -  (^  +  ^)Ad2  +  ^  A(K^.  V)2]p 


m    M        m  M 
m  M   ^    x^  W 

and  this  gives  the  result  for  j3  "because  again  all  the  source 
terms  are  differentiated  at  least  once  with  respect  to  either 
y  or  z. 

Finally,  equation  (5B)  may  be  written  in  the  form 
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The  terms  on  the  right  side  are  all  smooth  across  the 
planar  segment  because  we  have  shown  that  a  and  p  are 
smooth  and  differentiation  of  Vl  =  Y  with  respect  to  z  or  y 
removes  any  possible  discontinuity  of  y.      This  implies  that 
the  function  y        is  smooth,  which  is  possible  only  if  y   is 
smooth.   Therefore,  the  planar  segment  for  flow  past  a 
source  is  not  a  surface  of  singularity. 

C.    The  Solution  Near  the  Remainder  of  the  Characteristic  Cone, 

The  teims  in  the  solution  which  are  singular  on  the 
characteristic  cone  are  those  which  are  superpositions  of 
plane  wave  functions  with  real  wave  speeds.   These  terms 
form  the  function  J,  (equation  (-43)  )•   Triey  are  similar  to 
the  terms  which  comprise  the  plane  wave  representation  of 
the  fundamental  solution  of  the  Cauchy  problem  for  a  hyper- 
bolic equation.   The  behavior  of  the  fundamental  solution 
near  the  characteristic  cone  has  been  treated  in  the 
hyperbolic  case  [  3  ],  Our  discussion  will 

therefore  be  brief  and  will  serve  chiefly  to  show  the 
similarity  with  the  hyperbolic  case. 

We  consider  the  solution  in  the  neighborhood  of  a 
point  Pq  =  {xq/zq,    Yq/zq)   which  lies  in  the  interior  of  a 
smooth  segment  S  of  C. .   We  assume  S  does  not  intersect  any 
other  singularity  surface  and  that  the  curvature  of  C^  at 
"p^  is  not  zero.   As  noted  in  Section  5-F,  it  is  possible  to 
draw  two  more  tangents  to  S  from  points  sufficiently  close 
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to  one  side  of  S,    the  plus  side^  than  is  possible  for  points 
near  the  other  side^  the  minus  side. 

We  rewrite  J  in  terms  of  integrals  over  0  <_  9  <_  tt 
rather  than  0  _<  9  <_  2^,    using  equation  (45).   Since  p"^  e  S, 

we  have  g^CpQ^SQ)  =  g^  e^^o^^O^  "  °  ^°^  ^°"^°  ^0  ^  t^^^]- 
There  are  two  complex  conjugate  roots  of  g.(pj0)  near 
p  =  p^  and  9   =   9    ,for   p  on  the  minus  side.  These  merge  on  the 
real  axis  at  0  =  0„  when  p  =  p^^  and  give  two  real  singular 
points  for  p"  on  the  plus  side. 

We  first  consider  a  point  p"  on  the  plus  side  and 
near  "p^.      In  a  neighborhood  of  9^,  there  will  be  two  real 
singular  points,  a-,  and  cc^.      We  move  the  contour  from  the 
real  axis  to  the  contour  c',  as  in  Section  6-A.   The  contour 
c'  lies  on  the  real  axis  for  9   near  9        except  for  semicircles 
in  the  upper  half  plane  of  radius  e  about  9  =   a-,    and  9  ==   ap. 
For  p  sufficiently  close  to  p"^,  we  may  move  c'  to  a  contour 
having  one  semicircle  in  the  upper  half  plane  which  avoids 
the  points  9^,  a-,   and  a^,  instead  of  the  two  (smaller) 
semicircles.   Furthermore,  this  semicircle  may  be   so  chosen 
that  it  avoids  these  three  points  as  "p^  -^  p"-.   This  follows 
from  the  continuity  of  the  singular  points  in  p"  .   We  thus 


find  the  solution   and  all  its  derivatives  remain   finite 

■> 

0 


as  p  — >  p„  from  the  plus  side 


We  now  let  p   be  a  point  near  p„  on  the  minus  side. 
There  are  two  complex  conjugate  roots,  p,  and  f3p,  of 
g. (Pj9)  near  9  =  9  ,  which  merge  on  the  real  axis  at  9  =  9 
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for  p"^  ■=  p"^.   We  move  to  the  contour  e'.   The  contour  c' 
avoids  all  other  singular  points  and  the  double  point 
0  =  0  in  a  neighborhood  of  "p^,  except  for  p.,  and  p^*   ^^ 
are  using  that  S  intersects  no  other  singular  surface. 
The  only  singularity  of  the  solution  as  p  — >  p^  will  then 
come  from  the  integration  over  a  small  neighborhood^  N(9q), 
of  9      on  the  real  axis.   The  term  in  J  which  will  be  singular 
as  p  — >  Pq  is 

-^  ^a+p+Y-3      p    z  g^lprs)  log  (zgi(pr9)) 

'i  ^  ^  ^a"l  ayP-l~a7^  J^/^  ~~~'\{^7^\{^)  ' 

a+  p  +Y  =  6  or  J, 

where  we  have  applied  the  Laplacians  under  the  integral  sign. 

We  assume  a  +  p  +  Y  =  6  and  apply  the  explicit 
derivatives  under  the  integral  sign.   The  behavior  of  other 
derivatives  of  I.  may  be  found  exactly  as  below. 


We  obtain 


,  ..'Y-1      „    cos'='"^0  sin^'^e  aT"^(0) 
Lz±L. Hfi — 


We  expand  g,.  {'p^O)    about  {'Pq,^q)}    ignoring  terms 
of  third  order^  to  find 

g.(pre)  =  r(9Q).(p^-  p^)  +n^0(0o)  •  (p-  Po)(^-eo)-  |R(^-^o)+- 

with  R  =  ir^(0Q)  •  (p^-  Pq)  +  A.^g0(0Q). 
We  then  have 

-8^- 


"z = L +F 


where  the  integral  of  F  remains  bounded  on  N(9^)  as  p  — >  p  , 
Neglecting  terms  which  are  finite  as  p  — >  p^  we  may 


'0^ 


approximate   I.    by 


,       .T-1  cos«-l   e      sinf^-1   9      ^Z"^(%) 

(48)    I.    ^  ^-^-p- ^-^ 

2Tr^z  Q^(9q,    -A.(eo)) 

J 


N(eg)[^(V("-"o)+ V(«o)-  (p"-ro)(9-9o)-  iRCe-V^l"^ 


The  vector  n  (^q)  is  a  unit  normal  to  S  at  "p^^  ,  since 
C.  was  constructed  (Section  5-F)  as  the  envelope  of  the  lines 
m.  n   which  intersect  the  normal  speed  locus  and  are 
perpendicular  to  the  vector  from  the  origin  at  the  point 
of  intersection.   The  tangent  line  to  G.  at  p  =  pp,  then 
has  n  (9p,)  as  a  unit  normal,  and   |  n  (9_)«(p  -  Pp,)  j  is  the 


perpendicular  distance  from  p   to  p^^.   We  note  that  this 


■> 

0' 

distance  cannot  "be  zero  for  "p  y  p'^, since  p   is  near  p^,  and  lies 
on  the  minus  side  of  S. 

We  now  show  that  the  constant  | RJ  is  the  radius  of 
curvature  of  C^  at  J^.      We  have  S^{'Pq3^q)    =  Tq'    ^""(^o^'^i^  ^0^ 
=  0,  so  R  =  A.  gg(0Q)  +  ^^(.^q)'      The  curve  G.  is  the  envelope 
of  the  curves  g-{p,d)    =   0  with  respect  to  9.      This  gives  the 


-85- 


simultaneous  equations  S,-{p,Q)   =  g^;  a{p,S)    =    0,    which 
may  be  easily  solved  to  yield 

p^=  A.  (0)  n^{9)    +   A.^0  n^0(^)  .     0  <  9  <  TT  , 
as  an  equation  for  C.   We  then  have 

K  r(0)  = ^  -  -  ^-^ 

where  k  >  0  is  the  curvature  and  N  is  the  unit  normal  in 


the  minus  direction^  giving  the  result.   We  note  that 

'o)  •  (r-  To'. 


R  n  (9^)   •  (p  -  Pp.)  >  0,    independent  of  the  relative 


directions  of  n  (^q)  and  N  (^q)- 

We  may  now  evaluate  the  integral  in  equation  (48) 
and,  keeping  only  the  highest  order  singularity  as  p  — >  p   , 
we  find 

I  ^  Lzl)l~l     cos""^  Qq  sin^"^  Op   Aj   (0o)  27r  sgn(Rl^^  ^ 
^  2Tr2z2    Q^(9q,-A.(0q))  (-2REr(r-  To)?"^~ 
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Figure  1.   Section  of  the  characteristic  surface 
containing  A^^. 


Region  I 


Figure  2.   Intersection  of  the  section  of  the  characteristic 
surface  containing  A„  and  the  tangent  cones  from 
-■u„.   The  shaded  regiun  is  the  planar  segment. 
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fast 


Figure  3.   Intersection  of  the  section  of  the  character- 
istic surface  containing  A„  (solid  lines) 
and  the  normal  speed  locus  (dashed  lines). 
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